CHARGE RENORMALIZATION AND STATIC 
ELECTRON/POSITRON PAIR PRODUCTION FOR A 
NONLINEAR DIRAC MODEL WITH WEAK INTERACTIONS 

JULIEN SABIN 

Abstract. The Hartree-Fock approximation of Quantum Electrodynamics 
provides a rigorous framework for the description of relativistic electrons in 
external fields. This nonlinear model takes into account the infinitely many 
virtual electrons of Dirac's vacuum as well as the Coulomb interactions be- 
tween all the particles. The state of the system is an infinite-rank projection 
satisfying a nonlinear equation. In this paper, we construct solutions to this 
equation, in the regime of weak interactions (that is, small coupling constant 
f~ l , a), and strong external fields (that is, large atomic charge Z such that aZ := k 

O f stays fixed). In this regime, we are able to remove the ultraviolet cut-off A 

as soon as a log A stays fixed. As an application of this result, we compare 
the critical strength Kc(oi) of the external potential needed to produce an 
additional particle in the vacuum, when a = or « > 0. We prove that 
lima_>o K;c(a)/Kc (0) > 1, and we identify the limit exactly. Because of the 
dielectric behavior of Dirac's vacuum, static electron/positron pair production 
occurs in the interacting case for a stronger field that in the non-interacting 
case, which is a mere consequence of charge renormalization. 



Introduction 

The relativistic kinetic energy ol a quantum electron is described by the Dirac 
operator = ~ia • V -I- /3 acting on i^(R'^,C'*), where a. = (aj)j=i_2,3 and 



a, j' ^^-I'^'S. -Id. 

Ol\ / -i \ f I 



Here, the units are taken so that the reduced Planck's constant h, the speed of 
light c, and the mass of an electron me, are all set to 1. The spectrum of the 
Dirac operator is cr(Z?°) = (— oo,— 1] U [l,+oo). To explain why real electrons 
have a positive kinetic energy, Dirac suggested to postulate that, in nature, all the 
negative kinetic energy states of the Dirac operator are already occupied by virtual 
electrons. Thus, the relativistic vacuum is composed of infinitely many charged, 
virtual particles, usually called the Dirac sea. 

In an external electromagnetic field, the distribution of charge of the Dirac sea 
changes. This phenomenon is called vacuum polarization. In the same spirit, a 
negative kinetic energy virtual electron can turn into a positive kinetic energy real 
electron when given enough energy. This transformation results in a "hole" in the 
Dirac sea, which itself has a positive energy with respect to the fully filled Dirac sea. 
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This hole can thus be observed and is identified with the anti-particle of the electron, 
the positron. Hence, such a process is labeled electron/positron pair production. In 
this paper, we are interested in the case where this energy is brought to the vacuum 
by an external, classical electric field. This situation has been originally considered 
by Sauter[24], Heisenberg-Euler p2], and Schwinger [25]. Since it is non-negligible 
compared to the other sources of pair production only for huge electric fields, it has 
not yet been observed experimentally. However, the creation of such a strong field 
could be achieved in the near future, due to recent progress in laser physics [J]. 

There exists two approaches to the problem of electron/positron pair production 
under an external electric potential V. The first one consists in studying the vac- 
uum in its ground state, in the presence of V. One tries to determine if the vacuum 
then contains real particles or not, if V is strong enough. In other words, for large 
V one expects that it is energetically more favorable for the vacuum to contain 
real electrons or positrons, than to contain no real particle at all. If the vacuum 
contains real particles, we say that the external potential V produces pairs, or more 
precisely produces particles. Indeed, for this formulation of electron/positron pair 
production, called static pair production, it is not necessary that the number of 
electrons created is the same as the number of positrons created. We still refer to 
it as pair production though, because one can think of it as a localized particle 
together with an anti-particle at infinity. The second approach consists in study- 
ing the time evolution of the vacuum with no electron and no positron (the free 
vacuum), when the external field is compactly supported in time. One wants to 
determine if, for large times, some pairs still exist. In this case, notice that it is 
necessarily pairs that are produced since the total charge of the system is conserved. 
The second approach is thus called dynamical pair production. While dynamical 
pair production may be more relevant from a physical point of view, it is also much 
harder to study mathematically. 

Only a few mathematical results are known about electron/positron pair produc- 
tion. Indeed, one has to deal with infinitely many relativistic particles, interacting 
with an external field and possibly with each other. The correct theory to describe 
such systems is Quantum Electrodynamics, which has only been formulated in a 
perturbative fashion so far. Even the apparently simpler many-body theory of N 
relativistic electrons is not well-defined mathematically. For the time being, one 
has to make approximations to obtain results about such systems. 

The case where the interactions between the particles are neglected is one of 
these approximations. In this setting, everything reduces to the spectral theory of 
the one-body operator D'^ + V, for which the issue of dealing with an infinite number 
of particles does not exist anymore. The problem of static pair production has been 
studied by Klaus and Scharf |17] . The question of dynamical pair production has 
been tackled in the non-interacting case by Pickl and Diirr [20,, who proved that 
pairs were created for a certain class of potentials adiabatically evolving in time. 
This relies on a subtle analysis of the dispersive properties of the essential spectrum 
of the Dirac operator. 

When the interactions between the particles are taken into account, one can 
perform a mean-field approximation to again reduce the problem to the one-body 
space. The additional difficulty compared to the case without interactions is that 
the state of the system is not described by a wavefunction anymore, but rather by 
a one-body density matrix, an infinite-rank operator on the one-body space (usually 
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an orthogonal projection). The idea of performing a mean- field approximation of 
Quantum Electrodynamics goes back to Chaix and Iracane [2 . The first rigorous 
models formulating these ideas were developed by Hainzl and Siedentop p,4] . and 
by Bach, Barbaroux, Helffer, and Siedentop pj. Later on, another related model, 
which took into account vacuum polarization, was introduced by Hainzl, Lewin, 
and Sere in 2005 |10], and extended in a series of articles by Gravejat, Hainzl, 
Lewin, Sere and Solovej [HI [131 lEl E [71 [5] . In the present article, we will use the 
latter, called the Bogoliubov-Dirac-Fock model. In particular, the polarized vacuum 
is constructed in a rigorous fashion as a ground state of a certain nonlinear energy 
functional. Static electron/positron pair production has already been studied by 
the author in this framework in a previous work |23j . 

In the case without interactions, the one-body density matrix of the vacuum 
polarized by a potential V is given by 

If V is smooth enough, the eigenvalues of -I- kV are continuous with respect to 
the parameter k. Klaus and Scharf [17] proved (in the case F < 0) that a particle is 
created in the vacuum if one of these eigenvalues crosses when increasing k. If it 
crosses from the right to the left, it corresponds to the creation of an electron, while 
if it crosses from the left to the right, it corresponds to the creation of a positron. 
By the min-max principle of [5], in the case 1^ ^ eigenvalues can only decrease 
so only electrons are created. This picture was confirmed by Hainzl [9], from the 
point of view of the charge density of the system. When interactions are turned on, 
the one-body density matrix of the polarized vacuum satisfies a nonlinear equation 
of the form 

(0.1) P = X(-oo,o] (^° + ^^ + a^p), 

where Xp is an operator depending on P (see Equation ()1.2|) below) and a > is 
a coupling constant controlling the size of the non-linearity (or, equivalently, the 
size of the interactions). The number a can also be interpreted as the square of the 
charge of an electron, a = e^, in which case it should also appear in front of the 
electric potential V. Since we are interested in pair production in strong fields, it is 
also natural to consider a fixed potential V (interpreted as the shape of a nucleus, 
for instance), which strength is controlled by a parameter Z > (interpreted as 
the charge of the nucleus) . The equation satisfied by P is then 

(0.2) P = X(-co,o] + {c^Z)V + aXp) , 

As a consequence, studying ()0.1|) for V fixed and a > small is the same as studying 
()0.2|) for aZ fixed and a > small, which is a well-known regime in the physics 
literature [121 HH [SI [55] ■ In the sequel, we will thus write k — aZ. The existence of 
a solution to (j0.2p is non trivial. Actually, one can show that there is no projection 
on L^(M'^,C'') satisfying (|0.2I) . One way to circumvent this problem is to impose 
an ultraviolet cut-off A > 0, and to solve instead 

(0.3) P = X(-oo,o] (nA(i?° + ^iV + aXp)HA) , 

where Ha is the multiplication operator in Fourier space by 1_b(o,a)- We will show 
(Theorem [Ij that for a general class of V, Equation (|0.3p has a unique solution for 
K fixed if a is small enough, A is large enough, and a log A fixed. Furthermore, we 
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will show that this projection P satisfies 

P - X(-oo.o] (d° + ^ ■ . v) ^ 

V 1+ 3^^alogA J 

in the Hilbert-Schmidt topology, as a — > with a log A and k fixed. We interpret 
this result as a first order charge renormalization of our model, meaning that we 
can compute the screening induced by the polarized vacuum explicitly in terms of 
a and A. The term "first order" refers here to the fact that this is valid only in the 
limit a — >■ 0. Roughly speaking, if we put an external potential V in the vacuum, 
what we observe in return (in the limit a 0) is a screened potential V^cr, with 

(0.4) 14cr = tV. 

1 + l^alogA 

A similar formula was already obtained in [6], and it was valid for all a > 0. 
However, in [6|, the non-linearity Xp is simpler than the one we consider in this 
article (namely, it does not contain the so-called exchange term). Here, we obtain 
(I0.4p only in the limit a — > 0, but we manage to include the exchange term. 

Finally, we study static pair production for this model with interactions. Because 
of the nonlinearity Xp, it is not clear if the picture provided by Klaus and Scharf 
is still valid, that is if an eigenvalue of -D(a, kV^) := IIa{D'^ + kV + aXp)Ilfi^ crosses 
when k is increased, then a particle is created. Instead, we will compare two 
operators 

P+ = X(-oo.^+] (nA(i?° + kV + aXp^ )Ha) , 

and 

^- =X(-oo,M-] (HA(i?" + AcV^ + aXpjHA), 
where — 1 < < < < 1. If the operator 13° + kV has an eigenvalue 
X{kV) e (/z_, /Lt+) and a is small enough, P_ is interpreted as the polarized vacuum 
with a certain charge q, while P+ is the polarized vacuum with charge q + 1. We 
want to determine which one is the most advantageous energetically. To do so, we 
introduce the function 

F{K,a)^£{P+)-£{P-), 
where £ is the energy functional which will be defined in the next section. Our 
second result (Theorem [5]) provides an explicit formula valid for a small enough 
and a log A fixed 

F{K,,a) = A(KV;er) +Ca^o(a), 
where V^cr is given by (10. 4p . In particular, if X{kV) is positive for k < Kc and 
negative for k > Kc, then the polarized vacuum has charge q for k < Kc{a) and 
charge q + 1 for k > Kc(a): an electron is created; if X{kV) is negative for k < Kc 
and positive for k > Kc, then a positron is created. Furthermore, we have 

Kc(a) — ^ (^^^ ^"^°^^) ^ 

as a —7- with a log A fixed. Hence, due to the screening of the polarized vacuum, 
static pair production in the interacting case requires a stronger potential than in 
the non-interacting case. 

In principle, we can expand explicitly the function F to any order. For instance, 
we can compute the next order of F in a, that is \ima^oiF{K, a) — A(KT4cr))/a. We 
are interested in the sign of this quantity, which would determine the growth of F 
to the next order. Unfortunately, we were not able to determine its sign, due to the 
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complicated form of this term. It remains an interesting open problem to compute 
this sign. Finally, let us mention [27], where the same regime a — > 0, a log A fixed 
is studied without external field and for the vacuum of charge 1 , interpreted as an 
electron together with the Dirac sea. 

The paper is organized as follows. In Section 1, we introduce the model and 
present our main results. The proofs are given in Section 2. 

Acknowledgments. I am indebted to Mathieu Lewin for his advices and his 
help. I also acknowledge support from the ERC MNIQS-258023 and from the ANR 
"NoNAP" (ANR-IO-BLAN 0101) of the French ministry of research. 

1. Main Results 

In this section, we first explain the setting of the Bogoliubov-Dirac-Fock (BDF) 
model, and we state our results later. 

1.1. The Bogoliubov-Dirac-Fock model. The content of this section can be 
found in [lOllIIlITSlliaE]. Let L^{R^,C'^)- For any A > 0, we define 

^A:={/e-^, supp/c B(0,A)}, 
where the hat denotes the Fourier transform of /, defined for / integrable by 

and B{0, A) denotes the ball of center and radius A in C^. We call such a A a cut- 
off in Fourier space. Since the Dirac operator is a multiplication operator by 
p 1—^ Ck-p+P in Fourier space, it stabilizes f)A and its restriction to Sja as an operator 
is denoted by . Notice that D° is bounded and that ||-D5(|| = Vl + A^. We 
use the notation _ = X{-oo,o] {D%) (resp. P?_+ = X[o,+oo) for the negative 

(resp. positive) spectral projection of the Dirac operator For any bounded 
operator A over j^A and any e,e' = +, — , we define the block matrices A^g' :— 
P° gAP°_g,. We also define the associated spectral subspaces iOA,± ■— P\^±^a- For 
any separable Hilbert space H and any p > 0, we denote by GpCH) the Schatten 
class of all bounded operators A over H such that \\A\\q :— Tr \A\p < +oo, where 
1^1 — A* A. We then introduce the vector space 

<3i,po _ (^a) ■■= {Q e &2{^a), Q++, Q— e Gi{^a)} ■ 

It is a Banach space endowed with the norm 

\\Qhpo^_ ■■= \\Q++\\e, + \\Q—\\e, + \\Q+-\\e, + \\Q-+\\e,. 
For any Q di &i po , we define its generalized trace by 

Tro(Q) :=Tr(Q++ + g__). 

Any Q S &i,po has an integral kernel Q{-,-) £ L^(M.^ x R^,M4{C)), where 
7W4(C) denotes the set of all 4 x 4 complex matrices. Thanks to the cut-off in 
Fourier space, it has a unique smooth representative allowing to define the density 
of charge pq associated to Q, defined by Pq{x) — Trc4 Q[x,x) for all x £ 'M? . It 
satisfies pq G L?{S?) n C, where C denotes the Coulomb space, 

C:={/, / 4#dp<+oo} = H-i(M-^). 
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For any f,g^C, the Coulomb inner product between / and g is defined as 



iiR3 br 

If /, g are smooth enough (e.g. L^^^ fl C), then 



/ 4^dxd,. 



y 

In the BDF theory, a quantum state is described by a self-adjoint operator P on 
Sja such that < P < Ijj^ in the sense of operators. Such a P is called a one-body 
density matrix. For instance, the operator P^ _ represents the quantum state of 
the (free) Dirac sea. For any one-body density matrix P, we define its renormalized 
one-body density matrix by Q = P — P° _ . The set of all admissible quantum states 
is then 

/Ca :={Qe 6i,po _, Q = Q*, -Pa,- ^ Q ^ 1 - Pa,-}- 

For any Q such that — P^ — PA_,we have Q++ > 0, —Q ^ and 

the inequality 

so that Q G ©2 if Q++,Q— G ©i- The quantity Tr(Q++) (resp. Tr(g— )) 
represents the average number of electrons (resp. positrons) of the system. Hence, 
the convex set /Ca represents the set of all quantum state with a finite average 
number of electrons and positrons. The BDF energy of any quantum state Q G /Ca 
is 

(1.1) 

£r^{Q) = Tvo{DlQ) + a [ pqV+^D{pq,pq)-^ [ f M^dxdj/. 

Jrs ^ ^ Jm3 jr3 \x - y\ 

Here, V is the external electric potential and a > is the coupling constant. The 
first term of this energy represents the kinetic energy of the system. The second 
term is the interaction energy between the external potential and the density of 
charge pQ of the system. The third term, usually called direct term, is the Coulomb 
electrostatic energy of the density of charge pQ. The last term, well-known in 
Hartree-Fock theory, is called the exchange term. It is well-defined for Q G /Ca by 
the Hardy-Kato inequality 

—- — — da^dy ^ -Tr(|P>A|Q )■ 
R3 \x — y\ z 

the right term being finite since is bounded. Since P^ ±-Da = ±1-0°! -Fa ±' 
kinetic energy can be rewritten as 

TvoDIQ = Tt\DI\{Q++-Q__), 

which is finite for Q G /Ca- The direct term is also well-defined since pg G C. 
Finally, we assume that V is given by the Coulomb potential generated by a fixed 
external density of charge G C of strength Z > 0: 

SO that we can rewrite 



/ 



PqV = -D{pQ,Zv). 
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With this choice of V and by the Hardy-Kato inequality, wc find that £^^^ (Q) is 
well-defined for any Q e /Ca, and that it is bounded-below for ^ a ^ 4/7r since it 
satisfies 

It was proved in [11] that for any i' E C, for any Z > 0, and for any ^ a < 4/7r, 
the BDF energy has at least one mininiizcr on ICa- If = 0, then Q = is the 
unique global minimizer of the BDF energy. If 7^ 0, uniqueness does not hold 
in general. However, there always exists one minimizer Q such that Q + _ 
is an orthogonal projection. More precisely, this minimizer satisfies the following 
self-consistent equation 

n o^ ( Q = X(-oo.o]{Dq) - Pi-, 

where IIa is the multiplication by 1_b(o,a) in Fourier space, Vp := | • for any p, 
and Rq := for any Q. This minimizer is interpreted as the polarized vacuum 

in the presence of V,y. 

1.2. Static Pair Production. We now present the setting we use to study pair 
production. As we said in the introduction, we consider the vacuum polarized by 
a density of charge Zv. Moreover, we consider the limit of weak interactions, ie 
a — ?> and strong external field Z — >■ -l-oo, with aZ fixed. Hence, we define a new 
parameter, 

K := aZ. 

We make the following assumption on v G L^(R^) n C. We assume that there 
exists an eigenvalue A(kj/) of multiplicity 1 of the operator D'^'^ := Z)° — kVi, on 
that crosses at some value k = Kc G M. In mathematical terms, it means that 
there exists Kc € M and e > such that 

• For all KG [kc — s, Kc + £], the operator D'^'^ on Sj has only one eigenvalue 
X(Kh') in some fixed neighborhood of 0, of multiplicity 1; 

• For KG [kc — e, k^), > 0, and for k G (k^, + £]j A(ki/) < 0. 

The assumption of multiplicity 1 is made for simplification and it can be easily 
removed. For the sake of clarity, we will also absorb He in ^ by assuming in the 
following that 

Kc = 1. 

We could also assume that \{kv) changes sign the other way around, which only 
changes the physical interpretation: a positron is created instead of an electron. 
Of course, since v G nC, the eigenvalue \{kv) is a continuous function of 

K. Hence, there exist — l</i_<0</i+<l such that for all k G [1 — e, 1 + e], 
A(k) G (At-,M+) and n± ^ a{D'"'). 

The interacting vacuum polarized by Ziy is a solution to the equation 

^ ' \ Dq ^ HAiD-^ + a{Vp^ ~ Rg jjnA. 

Here, the operator Dq is seen an operator on S^a- The non-interacting vacuum in 
presence of the density of charge Zi/ is given by 
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where = X(-oo,o] is the negative spectral projection of the Dirac operator 
in ^. Notice that, with our notations, the non-interacting vacuum satisfies 

/ Qii„(«) = Qiin,-M = xi-ocM-]!^""^)-^- if «:<i; 

1 Qlin(«) = Qli„.+ («) = Xi-oo,^.^]iDn - P- if '«>1- 

Consider the non-hnear analogue of these equations, ie 

(1.4) Q±{K,a) = X(-oo,M±](nA(i?"'' + a{Vp^ - Rq^a) - P^._. 

Our first result shows the existence of solutions to this equation. By an easy adap- 
tation of [TDl Theorem 2], it can be shown that, for a small enough, Q^{K,a) 
(resp. a)) is a minimizer of the BDF energy in the domain of charge 

Tro((5iin,-(K)) — q (resp. Tro(Qiin,+ (K)) = Hence, Q+{k, a) has one electron 

more than (k, a), as in the linear case. We want to know which one of these two 
has the lowest BDF energy. To answer this question, we introduce the function 

(1.5) F{k, a) = (Q+(«, a)) - £^^^ {Q_{k, a)). 

If F changes sign while increasing k, it means that the vacuum gains one particle, as 
explained in the introduction. Our second result provides an expansion of F{k,q) 
as a — > with a log A fixed. 

1.3. Statement of the Results. Our first result concerns the existence of solu- 
tions to the equation (|1.4p . In [TU], it has been solved by a fixed-point argument 
when /i = 0, for a small enough, for A large enough, and for v G Lp'iE?) n C such 
that ||Ki^||L2nc is small enough. Here, we need to solve it for fj, ^ a{D'^'^) and under 
no smallness assumption on the parameter k (indeed, if k is too small, D"" cannot 
have an eigenvalue near 0, which we assume in our application to pair production). 
Hence, we have to go beyond the method of [l^. A solution to (|1.3p can be found 
by a minimization argument, as in [11 . With this technique, no assumption on 
the smallness of IlKt^llL^nc ^-iid of a (except a < 4/7r) is needed. By the same 
arguments, one can solve (jl.4p with general fi by minimizing the BDF energy in 
sectors of charge N, given by IC\{N) := {Q e ICa,Ti-q{Q) ~ N}. The existence of 
minimizcrs in these charge sectors has been proved in [12^, under the assumption of 
a binding condition, well-known in many-body theories. However, we do not know 
if these binding conditions are satisfied in the regime of parameters we consider. 
Even if we could prove this binding condition, the range of fj,s for which we can 
solve (|1.4p . ie the range of the map iV i— >■ /i, is not known. This range has been 
studied carefully in [6] , as well as the verification of the binding condition, in the 
case where the exchange term is neglected in the BDF energy. In the setting of 
our study, nothing is known about this issue. For all these reasons, we solve p.4p 
by a fixed-point argument, following the method of [10]. This approach has also 
the advantage of furnishing a priori estimates on the solution of ()1.4|) in adequate 
norms, which are useful in proving the renormalization formula (|0.4p and also in 
the proof of Theorem [2l 

Before stating our result, we have to introduce the Banach space where we will 
find solutions to (|1.4p . which was introduced in [10]. For any a; G M'^, we use the 
notation E{x) \/l + x'^. We define the space 
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where 

Q:-{Qe©2(iDA), ||Q|||:= / / E{p-q)^Eip + q)\Qip,q)\^dpdq<+^}. 

The vector space A" is a Banach space, endowed with the canonical norm 

\\{Q,p)\\x:=\\Q\\Q + \\p\\L-nc- 

We find solutions to (|1.4p by applying a fixed-point argument to the couple {Q,Pq)- 
That is why the Banach space X contains both the one-body density matrices Q 
as well as the densities of charge p. 

Theorem 1. Let v e L'^{M?) nC, p e (-1, 1), and s ^ such that 
(1.6) Vk e [1 -e,l + e], /i ^ (£>"'')■ 

Assume furthermore that 

log(2-)- |fc|)2|9(fc)|2d/c < +00. 



Then, there exists Lq — Lq{p, u, e) > 0, ao = ctoip, v, e) > 0, Aq = Ao(/i, v, e) > 0, 
and Rq = i?o(/z, J', e) > such that for all ^ a ^ Uq and A ^ Aq satisfying 
a log A ^ Lq{p, v, e), and for all k £ [1 — e, 1 -I- e], the equation 

(1.7) Q ^ X(-oo,t.]{^A{D° - kV, + a{Vp^ ^ Rq))^^) - Pl_ 

has a unique solution {Q,Pq) £ X such that 

(1-8) \\{Q + Gifl{KUren), PQ - BAKiyren)\\x ^ Ro- 

Here, we used the notations 

(1-9) t'rcn := / p t^, 

1 + uBa 

and for any p G L^(M'^) H C, Gi,o(p) £ Q is defined by 
Finally, the function Ba is given by 

(1.11) BA = -r — ^d2 = — logA- — + — ^ + Oa^oo(A-2). 

Remark 1 . The function Ba already appeared in the previous works [101 1111 |B] . It 
is logarithmically divergent in the cut-ofF A. This kind of divergence is well-known 
is Quantum Electrodynamics, and has been noticed in the early days of this theory 
by Dirac. The reason why we study the regime where a log A is bounded is to 
control this divergence. 

Remark 2. One can apply this theorem to £ = 0, in which case the condition on 
p, V is simply p ^ a{D'^). Here, we allow e > since we need estimates uniform in 
K for the proof of Theorem O 

Corollary 3 (First order charge renormalization) . Under the assumptions of The- 
orem]^ and if we denote by Q{K.v,a) the solution of (jl.7l) given by Theorem]^ we 
have as a ^ with a log A ^ Lq fixed: 



Q{KV,a) - (9lin(Ki^rcn) > 
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in Q, where 



Remark 4. This corollary does not imply the convergence of Q{Kv,a) towards 
X(-oo,/i](^'*'^""') — -P° in Q as a ^ 0, merely because this last operator does not 
belong to Q. However, it belongs to 62 by [12] • Since convergence in Q implies 
convergence in 62, we deduce that 

in ©2, as a — > with a log A ^ Lq fixed. 

Remark 5. By (jl.lip . in the limit a — )■ with a log A fixed, we have 

1 



l + :^alogA 



3tt 

which was the form used in the introduction. 

Thanks to Theorem [TJ we can build a solution to (ITT71) with v e X^(R^) n C 
satisfying the assumptions of Section [TT^ for ^J. = fJ,±, for any kg [1 — e, 1 + e], 
and for any a > small enough, A large enough such that a log A is small enough. 
Since the corresponding Q±(kz^, a) belongs to Q C &2{^a) and is a difference of 
orthogonal projections, we can use the fact that 

to infer that Q±(k, a) G JCa, hence the function F given by (II. Sp is well-defined. 
Our second result is the following: 

Theorem 2. Let v G -L^(M'^) H C satisfying the assumptions of Section \l.S\ with 
furthermore 

log(2+ |fc|)2|9(fc)|2d/c < +00. 



Then, for a small enough, for A large enough such that a log A is small enough, we 
have the following expansion for all k, £ [1 — e, 1 + e] .■ 



(1.12) F(At,a) = A(KI/ren) + C'a^0(a), 



where Oa^a {a) denotes a function which absolute value is ^ Ca for a small, where 
C only depends on /i, u, and e. 

Remark 6. Theorem [5] shows that for a small enough, the interacting vacuum for 
K > Kc{a) has one electron more than the interacting vacuum for k < Kc(a), with 
Kcia) ^Q-fO 1 + ai^A > 1 = Kc- In words, charge renormalization hinders pair 
production. 

Remark 7. One could look for the next order in the expansion (|1.12p . i.e. look for 
a number c such that 

F{k, a) — A(Kz^i.on) + ca + OQ_^.o(a). 

Of particular importance would be the determination of the sign of c, since it would 
give the next order influence of the interactions on pair production. Indeed, the first 
term in the expansion of F, X{Ki'rea) is purely an effect of charge renormalization, 
and not a direct consequence of the interactions between the particles (though 
charge renormalization is itself a consequence of the interactions). It can be shown 
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that c ~ ((y5, iJ(Kfrcn)i^), whcrc G is the eigenvector associated to A(Kz^rcn), 
and H{Ki>,:cn) is a self-adjoint operator on Sj. However, we could not determine if 
H^Ki'ron) had a sign, because of its very complicated form. For this reason, we do 
not provide the next order expansion of F in this paper. The details are provided 
in |22] . In principle, our method enables to compute the expansion of F to any 
order. 

Remark 8. When the exchange term is included, the operators P" and should 
in principle be replaced by effective operators V'^ and D'^ on S]\ satisfying a non- 
linear equation and which depend on both a and A, see |13j . Our strategy works 
the same with V?^ and T)'^, but the estimates are even more technical. Some details 
on how to adapt our method will be provided in |22j . In particular, in Corollary [3l 
the operator D'^ — nV^^^^ becomes 

(1.13) Z?eff = g'Moc ■ i-iV) + 5o(0)/? - T-^K, 

where 

.i(0) = l + f«logA, .0(0)^-4^, ^^-f^- 
Stt l-3^alogA 37rgi(0) 

In terms of mass and charge renormalization, the interpretation of this operator is 
not completely clear to us (see [18] and [13l Sect. 2.5] for a discussion). 

The rest of the article is devoted to the proofs of Theorem [T] and Theorem [2l 

2. Proofs 

2.1. Proof of Theorem [H The strategy for solving (II. 7p in TU' was to see the 

operator 

Dq := Dl + aU^iVp^-Zu - Rq)^k 

as a small perturbation of D\ if a is small enough, and then use the resolvent 
formula to expand 

X(-oo,m](^q) - X(-oo,/.](£'a) 

in powers of a, to find a functional to which they apply the fixed-point theorem 
of Banach-Picard. Here, we obviously cannot follow this strategy since aZ = k 
is fixed and not small. Instead, we have to find a new reference operator -Drcfj 
independent of Q, satisfying Dq = Dj-cf + 0{a), in an adequate sense. We choose 



(2.1) D,,, Dl - HaK,_Ha + aHAi?i^o('«f^rc„)n^ 



wher£Qj^i,o('«^^ion) '■= R[Gis){ni^-[cn)\ is the exchange term associated to Gi_o('*'^rcn)7 
which was defined in p.lOp . and t'rcn was defined in (jl.9p . 

Remark 9. A more obvious choice for Z^iof would have been Z?icf = F)\ — Ha Vkj^Ha. 
Unfortunately, for this choice we do not have Dq — D-^d — C'(a) since there remains 
diverging terms in the expression Ii]^{VpQ — i?Q )Ha. These terms diverge precisely 
as log A, and extracting them leads to the expression (|2.1[) . 



^Here and in the sequel, for any operator Q with integral kernel Q(x,y), R[Q] (or Rq) is the 
operator which kernel is given by '^!^^'y^ i and for any density of charge p, V[p] (or Vp) is the 
Coulomb potential generated by p, i.e. V [p] := p * [ ■ | 
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With this new reference operator, Equation (|1.7|) is equivalent to 

(2.2) Q = X(-oo,H (^ref + [VpQ+p^^i - i?Q+Q„f] Ha) - Pref, 

where 

-Prof := X(-oo,M](Aof), Qrcf ^'rcf— -P?,- +^1,0 ('«^'rcn) , Prcf PP„f-P^ _-BaK1^„ 

Under adequate conditions on v, we will show that ||(Qrct, Prct)||A' — Ca--i.o(l), so 
that we are able to solve ()2.2p by a fixed-point argument. First, we need to derive 
some properties of the operator Z3rcf- 

2.1.1. Spectral properties of the operator i'rcf- We first recall the following general 
result, which was partly contained in [10 . 

Lemma 10. Let p G L^(M^) n C and Q e Q. Then, we have 

(2.3) \\Vp\\L^ ^ 2^/^||p|U2^c, 
and there exists C > such that for all 77 G K, 

(2.4) \\Vp\D^ + ^,,l^\y^^-£-^^y\\,. 
For any 77 G M, we have 

(2.5) lli?Ql^" + *^r'lle. ^;^^||g||Q, 
where 



(2.6) Co ^ inf sup f E(2xy f 



R3 ^(2u)i+-|7/-a;p 
Finally, we have the estimates, 

(2.7) PqI^VIU^^l^ = WRqWh^^l-^ ^ 2||Q||e„ 

(2.8) WRqIDT'^'Wl^^l^ = \\Rq\\hW2^l2 ^ |||Q||q. 

Proof. Inequality (|2.3p was proved in [Tni Lemma 3]. Inequality (|2.4p is a conse- 
quence of the Kato-Seiler-Simon inequality 

(2.9) Vp ^ 2, V/,g G Lf(M3), ||/(^)g(_jV)||e, (27r)-3/P||/||ip||.g||ip, 
of the Sobolev inequality HVpHie ^ C||VV^||i2 = C||/9||c, and of 



1 



\a ■ p + 13 + ir]\ 



6 



dp C 



Finally, inequality (|2.5I) is proved using the fact that for all p,q ^ 

E{p + q)E{rj) 



{E{pf + ri-^)E{p ^ q) 



^ 2, 



hence 



RQ{p,q)? 
e{pY + 
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where the second mequality was proved in [10, Lemma 8]. To prove (|2.7p . let 
(p £ . Then, by the Holder inequality we have for all x G R"^, 

\Rq^{x)\'^[ \Q{x,y)\'dy [ P^dy. 



\x - y|2 

By the Hardy inequality, one deduces 

\\RQ^\\h^m\%jy^\\h^nQ\\%j\\D>\\h- 

The proof of p.Sp relies on the same computation for (f e H^^^: 

Jr^ \x-y\ Jr3 \x-y\ 
so that by the Hardy-Kato inequality applied to both terms. 



□ 



Lemma 11 (Properties of Gi.o(p)). Let p e L^(R^) nC. Then, for any A > 1, the 
operator Gi,o(p) defined by p.lOp satisfies for some C > 0; 



(2.10) I|Gi^o(p)||q ^CVbiAllplU^nc, 

(2.11) ||Gi,o(p)||e. <4||p||c. 

Proof. Inequality p.lOp was proved in [101 Lemma 11]. It relies of the explicit form 
of the integral kernel of Gi.o(p) in Fourier space, obtained by the residuum formula: 

(2.12) Vp,<7eM3^ G^p){p,q)^-^^l^Vp{p~q)M{p,q), 
where the function M is defined by 

It was proved in [lOl Lemma 12] that 
(2.14) \M{p,q)\' ^16 



E{p + q)^' 
hence 

l|Gi.o(p)|||, ^ - / \mv,ik)\'dk[ ^ = l6||p||2, 

where we used that /jj,3 E{ey^ d^ = tt^. □ 

Combining Lemmas [TUl and [TTl we infer that for any a > and A > 1, the 
operator D^cf is bounded and self-adjoint on The following result ensures 

that Drcf has a gap around p for a small enough and A large enough. For any 
uj £ i^(]R'^) n C, we introduce the auxiliary operator on Sjj^, 

(2.15) Aef(^) -Da - HaVLHa + aHAi?i,o(w)HA, 

so that Drcf = Drcfini^Tcn)- 
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Lemma 12. Let v e L^(M^) n C and ^ e (—1, 1). Let e ^ such that for all s e 
[1 ~ s , I + e], ^ ^ (7sj [D"'^). Then for all L > 0, there exists ao ~ ao(/i, i', e, L) > 0, 
Ao = Ao(/i, z^, £, L) > 1 and = Co(a'i > smc/i that for all ^ a ^ ao, 
for all A ^ Ao satisfying a\/TogX ^ anrf /or aZ/ s G [1 — £, 1 + e], w;e /lave 
d{fJ.,(J^j^{D,c{{si^))) ^ Co- 

Proof. The proof is done by contradiction. Assume there exists q;„ — > 0, A„ -> +oo 
satisfying a„VlogA„ ^ L, s„ g [1 — e,H-e] such that ci(/L(, cr^^^ (-DrefCsn^^)) 
Up to extracting a subsequence, we can assume s„ — >■ s S [1 — e, 1 + e]. Since 
(i(^,CTi3A^ (A-cf(srii^))) ^ V") there exists (/3„ e io„ := i^A,, with ||</Jn||L2 = 1 and 
\\{Dn — ^J■)^Pn\\L^ ^ 2/n, with £)„ :— I?rof(s„i^). For any 77 G M and ri G N we have 

= (1 + n„(a„i?„(s„J/) - K„<.)(-D° + ivr^){D° + iri)(pn - irjipn, 
where we used the notations n„ :— IIa,, and i?„(s„J^) := i?i,o(sni')- Now for 77 large 
enough, we have \\Vs„^{D° + ir])-^\\ < 1 by (EH), and \\anR^^^{uj){D° + i7j)-^\\ < 1 
by (|2.5[) and the fact that \/log A„ ^ L. Then, we have for such a 77, 

{D" + i7;)(^„ = /, , Tj t — B~f — r~~f7 — u7;o~r~"i^(^"'^" + *^'^")' 
(1 + n„(a„i?„(s„i/) - v;„^)(D" + 777) ^) 

showing that {ipn)n is bounded in i7^(R'^). Hence, up to a subsequence we can 
assume (pn ^ ^ i'n H^(M.^). This imphes that D^ipn D^ip and /j,(p„ flip weakly 
in L^(R^). Since Vsv{D'^)~^ is a compact operator we also have Vs„,j(pn Vs„(p 
strongly in L^(R^) and hence UnVs^uifin -> Vs„(p strongly in L'^{R^). Using (|2.7p 
and p. lip , we infer that 

i')\D I ||i2_j.2,2 ^ 2a„||Gi^o('Sri'^)||e2 ^ 8a„(l + e)||i^||c ^ 0. 

We thus have n„a„i?„(s„i^)<y9„ converging strongly to in L^. Hence, Dnfn ^ 
{D^ — Vsu)^^ weakly in and since ||(£>„ — /i)(/9„||i2 < 2/?t,, we have {D^ — Vsu)p^ — 
flip. But fi ^ a^-j{D^^), hence we must have = 0. This implies that UnVsni^Vn ~^ 
strongly in L^(M^). We deduce (1?° — fi)(pn strongly in L^, which is a 
contradiction by the Weyl criterion for the essential spectrum since |j(/3„||x,2 = 1, 
(fin ^ weakly in L^ and /i G (—1, 1). □ 

Lemma 13. Let lo G L'^{R^) n C and ^ G (—1, 1) such that ^ Dt-c{{uj). Then for 
all a ^ such that 8a||a;||c < 1 — |/i|, the operator \D'^\{D^cti^) ~ l^-)^^ (extended 
by on S^j^) is hounded on L'^{R^ ^ C**), and we have 

1 



where 



-Dref(w) - 



l-\^l\J 1-ImIV d(M,^(A-cfM)) 

Proof. By the resolvent identity, we have 

1 1 1 1 



n f ^ no ' Tfi nA(K; -ai?i,o(t^))nA- , , 

Dref(wj-/i T>A^t^ l-^ Urcf(UJ) - ^ 

Using |||i?0|(DX-/i)-i||i2^i2 sc: (I-ImI)"^ and \\R,,o{uj)\D"\-Hl-^l- ^ S\\uj\\c 
by pT7)) and (pTTTj) . we obtain 



CHARGE RENORMALIZATION AND PAIR PRODUCTION 



15 



□ 

Lemma [13 ensures that we can use Cauchy's formula to P^-of since D^^f has a gap 
around /i. The next result is the analogue of [10, Lemma 9]. It ensures that the 
operator Dq also has a gap around fi. For any Q ^ Q and oj, p' E L'^(M.^) DC we 
define the following operator on 

(2.16) D^,p,,.Q = Acf(w) + (Vp, - Rq) Ha. 

With this notation, we have Dq = pQ+p,.ef,Q+Qrcf 

Lemma 14. Let oj e L'^{R^) DC, p. e {-1,1), a ^ 0, A > 0, (Q, p') e X be such 
that the following conditions are satisfied: 

(1) M ^ ^(A-cf(^)); 

(2) 8a||^||c < 1- ImI; 

(3) a6||(Q,p')IU < 1, 

where, using the notations of Lemma \13l 

d[p,a(D^cf[i^))) 
Then, we have the operator inequality 

\D^.p',Q-f^\ ^ (l-a6||(Q,p')IU)|^rcfM-Ml- 
In particular, p ^ o'(£'lj.p',q)- 

Remark 15. The conditions listed in the above lemma give a non trivial relation 
between a and p (recall that D^cfiuj) depends on a). Fortunately, these conditions 
are satisfied ii p G (—1, 1), p ^ a{D^ — VL,) and a is small enough by Lemma [T^ 

Proof. We adapt the proof of W, Lemma 9] to our case. First of all, since 
crcss(A'of(t^)) = [-E{A),-1] U [1,E{A)] and p € (-1,1), p i ^(L'rcf (t^)), we have 
d(/x, cr(-Drcf('^))) > 0. Let u e i^A- By Lemma[TUl we have 

\^kVM\l- WAlAAl- ^ ■f'^''^"'^;^' ll|A-ef(^)-A^|^|U^, 

where we used the operator inequality |Z)rcf('^) ~ lA^ d{p,a{D^cf{A))- For the 
exchange term, we write 



\\UaRqu\\l2 ^ WRglD'^AWe 
By Lemma [T^ this implies 



|-Drcf('^) - P 
TIaRquWl: ^ Co\\Q\\Q0ip,Uj)\\\D,AA~f^ML2 



lll-Drof(w) - p\u\\l^. 



showing that 



|nA(v - i?Q)nA| (-, — . ... + 9{p,a) ||(g, p')IUlA-cf(^) - /^|. 

\d{p,a{Drci{uj))) J 



□ 
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2.1.2. Resolvent expansion. Now that we know sufRcient conditions on our param- 
eters to prove that both Z^rof and Dq have a gap around fi, we can use Cauchy 
formula, and then resolvent expansion (at least formally for the moment), to see 
that (|2.2p is equivalent to 

277 V-Dref-M + i?? Dq - fj, + ITj J 

which leads to 



E 



„>;^ 27r 7r Drci - fi + ir] 



1 

■ Drcl - 11 + 111 



We write it as 

(2.17) Q = E [Q + Qret, + Prof), 



where for any g e Q, w, p' e L^(M'^) n C, 
(2.18) 

f-l^n+l /■ 1 / 1 , 

= W — / TT^^ — nA(V - ^Q)nA7r^^ r- d,?. 



27r Jr L'rcf(w) - fi + iri \ -Drof('^) - M + 

We will prove in particular that (|2.17p makes sense by showing that the series 
converges in Q, for a small enough. Equation (|2.17p implies that 

(2.19) PQ = E «>[Gr™(Q + Qrcf,PQ + Prof)]. 

We use the same trick as in 10 , by writing 

p[G^iQ,p')]^p[GMp')] + JiiQ,p'), 
where we recall that Gi_o(p') was defined in (ll.lOp . It was proved in [10] that 

Vk e M^ J- [p[Gi,o(p')]] (k) = -BA{k)p\k), 

where J^[f] := /. The function B\{k) is computed in [5J Eq.(88)]. It satisfies 
B\{k) = for |fc| > 2A while for |fc| ^ 2A we have the formula 



1 /^-'"-D z'~zV3 



2ttJo VTT7^-\k\z/2 



with 



%ATA^-v/l + (A-r)2 
ZA[r) = , Vr ^ 0. 
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Notice that Ba = Ba{0). Defining J^[Q,p') := p[G';'-^{Q, p')] for n ^ 2, we see that 
(|2.19p is equivalent to 

V/c e M^ F[pQ + prcf](fc) = . ^ n /L^ ^f(fc) 

1 + aB\{k) 

Hence, the couple {Q,pq + Prcf) satisfies 

{Q,PQ+Pref) = $(g,PQ+Prcf) = (<3 + Qrcf , PQ + Prcf ) , *r"° (Q + Qrof , PQ + Prct) ) , 

where 

1 a" 

We prove that $ stabilizes a certain ball in X and that it is contractant on this ball. 
By the Banach-Picard theorem, it proves Theorem [TJ To do so, we need estimates 
on the function 

2.1.3. Estimates on $. We begin by a useful result about commutators. 
Lemma 16. Let p G L'^{^^)C\C and Q G Q. Then we have the following estimates: 



(2.20) V0<C<1, Vt>C + ^, \\[\D°\^,Vp]\D°\~%^^C\\p\\c, 

(2.21) V0<C<1, Vt>C, \\[\D°\'^,RQ]\DT'\\e,^C\\Q\\Q. 
Proof. We begin by noticing that, if T is a self-adjoint operator, we have 

\\[\Dr,T]\DT'\\B. ^ V2\\[\D'\i,T>]\D°n\e., 
where T>{p,q) — T{p,q)x\p\>\q\ for all p, q G M"^. Indeed, 



>\p\>\q\ ""IpKkl 

and using \T{p,q)\ = \T{q,p)\ we find 

ff \E{pi^m^^ ff |g(p)c-g(g)C|2 

— \T{p,q)\ dpdq. 

\p\>\q\ E(q) 

Hence, to estimate 1 1 [ | | , T] 1 1?° | 1 1 ©2 , it is sufficient to estimate \\[\D°\'^ ,T>]\D°\-^\\e2 
or ||[|DO|f,T<]|DO|-*|!e,. Now let A = [IZJOI?, 1/>]|D0|-*. For any \p\ > \q\, we 
have by the mean value theorem 

\E{p)< - E{q)<\ ^ CE{q)^-^\E{p) - E{q)\ < ^E{q)^-^\p - q\. 
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Hence, 



Pill. 



\p~q\^\Vp{p~q)\' 



dpdq 



^C\\p\\l 



since 2{t — C, + 1) > 3. To prove ()2.2ip . the same strategy only works in the case 
t ^ 1/2 ^ C (notice here that the equality t ~ ( ~ 1/2 is authorized). In the 
general case, first assume < C < 1/2. Then, using the formula 



VO < C < 1, > 0, x'^ 

we find 

{\Dr,RQ 

We then estimate 



sin(7rC) 



X + s / s 



ds 



sin(7r<^) 



(2.22) \\[\D'^\'^,Rq]\D'^ 



\DO\+s 
sin(7r^) 



[\D\Rq] 



ds. 



/•OO 


\D"\i 


Jo 


\D"\+s 



x\\\DTi[\D%RQ]\D'^\-i\\e, 
Setting B = {D'^l-^lD^l, Rq]\D"\^^ , one has 



X 

0|i-t 



\D^\ + s 



ds. 



\\Bf 



62 



E{p)E{q) 



dpdq 



^ C 



I 



E{p-qf 
K3xR3 E{p + q) 



\RQip-q)\^dpdq 



< CWQWI, 



where the last inequality was proved in [TOl Lemma 8]. Now assuming t > ( > 
we have r := ^ — i ^ 1, hence using 



(2.23) 
we deduce 



Vr s; 1, Vs ^ sup 



1 



x^l X + S (1 + s)l-max(r,0) ' 

s'> ds 



6; 



s; c 



^ (1 -|_ s)min(t+1.3/2) 



IIQIIq, 



the right side being finite since t > ( and C, < 1/2. Ifl/2^C<1, we write 
\D°\<: = |i:)0|C'|i:)0|i/2 ^ith O < C' < 1/2 to infer 



0|C _ 



sin(7rC') f°° \D°\^/^ ds 



iDOj + ssi-C 



Now for all s ^ we have 

[p^'^a] = p^(^[|^°|^^^^aH[l^T/^li:>°|i^al^°lDp^ 
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If Bi = |D0|-i/2[|7^0|3/2^^<]|£,0|-i^ then 

JJ\p\<\q\ ^[P + Q) 

If B2 = \D°\-^[\D°\^/^, |D0|i?<|D0|]|i?0|-3/2, then 

,2 _ /r |i?b)l/2-i?(g)l/2|2^(p)2^(g)2 _ 



We deduce that if B3 = [jD"!?, i?<]|i:)"|~*, we have, using (P?^ again 



llSslle. ^ CWQWq I 

JQ 



'\DO\ + s""\DO\+s |DO| + s"" \DO\+s "J s^-<' 

< CWQWq I (1 + 5)min(i+l/2,3/2) + (1 + 5)min(t-l/2,l) ) '^^^ 

the right side being finite since t > C' + 1/2 = C and 3/2 — C'>1- Q 

The next proposition is the analogue of jTUl Proposition 10]. It is the main tool 
to prove the contractivity of the map 

Proposition 17. Let uj G L^{M!^) DC, ^i e (-1, 1), a ^ 0, A > 1, (Q, p') e X be 

such that the following conditions are satisfied: 

(1) p i a(A-cf(^)); 

(2) p i <j{D^,p',q); 

(3) 8a||^||c<l-|Ml- 

Then, we have the following estimates: 



II /), ^2iQ, P'Mx ^ llPref IL.nc + c^V^^ IKQ, P')WX 

+ ^(aS||(Q,p')ll;,)", 



ri-l 

X 



(2.24) ||(($-)'(Q,p'),(<i>2)'(Q,p'))K«Vl^S + a5]n(aS||(Q,p')l 

n^2 

(2.25) 5 C(a.)(1 + a^l^)'" (l + lkll"nc) ^(m,^)', 

where C{p) > is some constant only depending on fi, and 9{fi,uj) is defined in 
Lemma \13l 

Proof. We start by estimating = X)ri>i '^"^n norm || ■ || q. The difficulty relies 
in the fact that the definition of the norm || • ||q uses the Fourier representation 
of the kernel of the operator considered. In [TD], it was not an issue since all the 
studied operators, for instance the resolvent (D^ +ir])~^, have an explicit kernel in 
Fourier space. This is not the case in our setting: the kernel in Fourier space of the 
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resolvent {Drcf{^^) ~ ^J■ + iv)^^ is not explicit. It is thus harder to estimate the Q- 
norms of these operators. To circumvent this difficulty, we will expand the resolvent 
(£'rof('^) ~ M + i'n)~^ in terms of the resolvent {D\ — /i + iri)~^ and of lu, in each 
separately. The terms only containing resolvents of the form {D\ — /i + ir])~^ can 
be estimated via the results of [10], while the terms still containing resolvents of 
the form (Z?iof(w) — M + j??)"^ decay faster and can be estimated using inequalities 
of the form ||Q||q < || |£'°|^Q||e2 • Indeed, the norm HlD^I^ • ||e2 is a rough upper 
bound for the Q-norm, but has the advantage of being computable without using 
the kernel in Fourier space of the operator. We thus expand the resolvent as 

^ ^' n f / = ^o('?) + ^i^'?)' 

Drcf[UJ) - fi + n] 



with 



(2.26) A,irj)^{ . \ . Yiu;,a^ ^ 



Let n ^ 1 and write 



Y{Lv,a) :=nA(K;-ai?i,o(a;))nA. 

(_'l\n+l r " , , 



27r ^ 

2:e{oa}"+i 

where we used the notation X = — Rq. The term G„^(o,...,o) was already 
estimated in [ini Proposition 10], 

(2.27) ||G„,(o....,o)||q ^n{l + (av/b^)^)"+^(l + ||^|li.nc)"+'ll(Q,P')r, 

where k„ = ^(/i)" for n ^ 2 and ki = C(/i)-\/log A, for some number C(/Lt) only 
depending on ^. Actually, the estimate (|2.27p was obtained for /i = in fnT. The 
way to obtain it for any ^ is to replace Equation (58) in [10] by 
(2.28) 



^(^^(p) - pY + - P? + E{p + q)E{ri) 

for some C'{^). It remains to treat G„_a; with a; ^ (0, . . . , 0). Using that for any 
p, g e R'^ one has 

E{p - qfE{p + q)^ E{pf + 2E{pfE{q) + 2E{p)E(qf + E{qf , 

we deduce that 
(2.29) 

IIQIIq < + 2\\\D°\Q\D°\m^ + 2\\\D^\^Q\D^\\\l^ + \\Q\D^\i\\l,, 

for any Q € Q such that the right side is well defined. To estimate ||G„.:i;||g, we 
will estimate the right side of (|2.29p with Q = Gn,x- 
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Lemma 18 (Estimates on v4o(?7), ^1(77)). Let 77 e M. We have the following 
estimates. 



(2.30) 



\A,m + II 1^1(77) II < c{fi){i + {aVW^r)M%^c(^i^,,L,) 



For all Q ^ a < ^ , we have 



(2.31) 



and for any B G ©2 such that \D^\i^" B G ©2; we have 
(2.32) 

\D'\i+'^A,{^)XB ^C{^l){l+a^/^t{l + \\uJ\\l.^c)m,P)\U |i?"|^+^i? 
©2 

We also have 



(2.33) 



Ai{il)\D'\i 



62 



C{iJi){l + a,/\^y\\uj\\l.^cOM 



and for any B €&2 such that B\D°\2 £ 62, we have 
(2.34) 

BXAo{v)\D"f^ ^ <C(/x)(l+aytoiA)Ml + ll^lli^nc)ll(Q>P')lk ^^l^^l^ 
©2 

Finally, we have 



(2.35) 



|l^"l^i(^)| 



©6/5 



and for any B such that \D'^\B G ©6/5' have 
(2.36) 



|i?°|^o(r])^S||e,,. =^C'(A*)(l + «Vl^)'(l + ll^lli^nc)ll(Q,P')IU \D'\B 



©6/5 



Proof of LemmaUE To prove ([230| . just use || |i:)°|(i:)^ - A* + ivy^h'^^L^ ^ 
C{^i), \\Y{u,a)\D^\-^/^L-^L- < C{l + aVW)\\uj\\L-nc and || |I?"|(A-cf(c^) - A' + 
«^)"^IIl2^l2 s; e{n,uj). To prove (lOTD . we write |i:'°|i+'^Ai(77) = ^1(77) + B2{v) 
with 



\D"\ 



D\ — ^ + ir] D\ ~ ^ + ir] 



1 



\D°\ 



0|i+(T 



— li + irj 
1 



1 



Drcf (w) - fi + irj 
4 



^A - M + 



r(w,a; 



1 



Aof (w) - fi + iri 



Now Si (77) = -811(77) + -612(77) with 
l^°l 



i?ii(^) 



- fi + if] 

X (Y{LO,a) 



-^A — f- + 



{Dl - + 777) 



Dl-fl + i7j' 

1 

-Drcf (w) - H + irj' 
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Bi2{r]) 



\D'\ 



D\ - fl + IT] 

|Z?0|t 



1 



for some i +a <i < 1. Using |L>0-Ai+ir/|-i ^ C(Ai), \\Y{ijj,a)\D°\-'2\\^2^^2 
C(/z)(l + aVbp:)||w|U2nc by ^ and (HH), - /i + z77)(A-cf(c^) - + 

^ C{fi)9{pi,ijj), i + (T < 1, and the two inequalities 



(2.37) 

obtained by ((2^ 

(2.38) . 

iJ^' — fi + ir] 

obtained by and (|2.10p . we get 



i?l,o(t^)- 



1 



£;(r/)i/2 ' 



Cill)y/\0gA\\uj\\L2r,c 



(2.39) 
hence 



1 



< 



6e 



\Biime, ^ 



C(/i)(l + aVbp:)'||w|li2nc^(M,'^) 



£;(r/)3/2 

The term Bi2{rj) is treated the same way, except that we use the commutator 
estimates of Lemma [Tni to get 



[|i?°r,y(c.,a)] |i^"r*||e. ^C(/i)(l + aVl^)ll^llL^nc. 



We thus find 



\Bi2{v)\\&, ^ 



C(/i)(l + aybp:)'l|c^l!i2nc^(M,^) 



We next set i?2(?7) = B2i{r]) + B22{'n) with 
l^°l 



B2i{ti) 
and 

B22{V) 



D\ - pL + it] 



J U A li n~ I'll 



1 



-D° — f-L + ir] 

|£,0|u 



,0|i+<T 



-^A - fJ- + "ill 



',F(u,q) 



M + ir] 



for some 1 + tr < s < | and s — ^ < u < 1. Such a couple (s, m) exists since a < ^. 
Using the estimates as before, we obtain 

\\B2{v)\\e, ^ 



C{]i){l+a^T^nu;\\l,^c(^{]i,u;) 
E{r]Y/^ 
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We thus have proved ()2.31|) . Let us now prove (|2.32[) . Hence, let B E &2 such that 
^\i+"'B E &2- RecaU that Ao{'ri) is the sum of five terms. We will treat the first 



\_u"\ ^ ' D E '02- fiecaii inai 
two independently. We thus write 

The first term to estimate is 



Hi := 



D\ — ^ + iri D\ — ^ + ir] D\ — fj. + ii] 



By Lemma [T6l we know that 
(2.40) \\X\D°'-^ 



\D°\i+'',x] ID' 



'\D"\^+'',X 



B. 



^C\\{Q,p')\ 



X, 



|0|-|-c 



Hence, 

Next, we treat 

H2 



^C\\iQ,p')\\x. 



\m\e,^C{f,mQ,p')\\;,\\\D^\i+-B 



l+cr 



Y(uj,a)—^ — XB 



+- 



\D°\ 



— H21 + -ff22- 

The first term is easily estimated by 



J D^ - 1.1 + ti] 



H2i\\e, i^Cifi)\\Y{uj,a)\D\ „„ ,„ „ 

C(/i)(l + av/biA)||c.||i2nc||(g,p')lk 



s^C(/i)(l + av/logA)||c.||i2nc||(g, 
while we notice that the second term can be estimated by 



ID^I^+^B 



62 



C(M)(l + aybiA)lhllL^ncll(Q,p')IU 



|D"|i+'^D 



62 



We conclude that 

\\H2\\e,^C{p){l+a./]^)\\u;\\L2nc\\{Q,p')\\x 
The last term to control is 



if3 := 
We use 



^A + fc^o 



Dj( - p + IT] 

that ||i/3||e2 < ||^^3i||L2^L2|l^32||e2, 



with 



|D"|i 



— t: I Y(uj,a) — n 
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On the one hand we have 



— IJ' + iV 



XB. 



||ff32||e. + av^)' (1 + ll^lli^nc) ll(Q,P')IU \\D'\i+'^ B 

while on the other hand 



H?,\ — 



\D°\ \D°\^+'' 



1 



- + 
^A — M + 



|D°|^+",r(w,a)l ID"!"* 



^A - + 

,0,-. I^°l 

- + i?7 



—Y{uj,a)\D°\-'^ 



|0|t-i 



- M + 



'a - M + 



for sonic 1 + it < i < | and < — i < s < 1 , so that writing 



ID"! |Z)°|H'^ 



1 



^A - M + 



.0,-1 W,. „M^O,-i 



- + i77 



lg°l 

£)0 - /i + 



r(tj,a)|i:'°r3 



^A - + «7? 
I^°l 



-D^ - + 



^A - M + 



r{- fi + ir] 



for some ^ + cr < u < 1, one has 



Hence, 



WHsWe. ^ C(m)(1 + av/biA)^ (1 + Mhnc) \\{Q,P')\\x \D"\^+''B 



which finishes the proof of p.32p . Let us thus pass to (j2.33p . The proof is the 
as p.3ip . with a additional commutator between jD"! 2 and (Drcf(w) — n + ii]) 
Indeed, we have 

Mv)\D'f^ = f „o \ - Yic^a))' \D'\^ / \D'\ 



same 
1 



1 



5 r 



We have by the same method as before 



D,-cf{uj) - fi + iri 



-AD 



,0|i 



Now notice that for any C G K, 
1 



Dref(w) - fJ. + ilJ 



lOiC 



Drefiuj) - H + if] 



Dref (w) - n + iri' 
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SO that T2 = T21 + T22 with 



-Yiu,a) \D'\ 



1 



Drcfiuj) - fi + irj 
X \DX' \\D"\KYiuj,a) 



\D°\x 



1 



1 



1 



-Drcf (1:^) - H + ir] 
\D"\\DT' \\D"\KY{u,a) 



-Drcf (1^) - H + irj 



Z^rof (w) - fl + il] 



\D'\ 



Drcfi^^) - fl + il] 

for some 1 < i < |. Hence by the methods introduced before, we deduce that 

leading to (|2.33p . The proof of ()2.34p is the same as the proof of p.32|l . Next, we 
prove (|2.35p . It is because of this term that we have to expand {Di-c{{io) — ji + irj)^^ 
to the 5th order, as we can see from the folfowing estimate: 

5 



|||i?°|^i(ry)||ee,, ^ 



{Dl-fi + 11]) 



Y{u:,a) 



1 



1 



I- fi + lT] 



&6 



C(M)(l + «VEp:)'ll^lli2nc^(M,c^) 



£;(?7)5/2 

which is exactly ()2.35p . Finally, to prove ()2.36p . we simply notice that 

\\\D'\Ao{v)XB\\^^^^ ^ \\\D'\Aom\X\D'n\\\\D°\B\\e,^, 



^ C(M)(l+avAogA)4(l + ||c.||i.nc)ll(Q>P')IUIl|i?°|i3|lee/., 
by inequalities (|2.30p and (|2.40p . This ends the proof of Lemma [THl □ 

Let us now estimate |||£'°|^G„,2;||62- First, assume that xq — 1. We have 

« n 

(2.41) \\\D'f-Gn.46.^ / \\\D'f^A^me.l[\\XA^^{i])U.^L.di]. 

We deduce from ([OOl) that for all 1 ^ j 71, 

\\XA^^(t])\\l2^l^ <: \\Vp>A^^{t])\\l2^l2 + \\RqA^^{ij)\\l2^l^ 

^ \\Vp'\\l-^l4AccAv)\\l-^l- + \\RQ\D°\~'\\eJ\D''\A,^{7])\\L-^L^ 

^ C(m)(1 + av/foiA)^(l + Ml2^cnfi,u;)\\iQ,p')\U. 

Hence, for all 77 € K, 
(2.42) 



'[[\\XA,^ifj)\\L2^L2 ^ (c(M)(l + av/fo^)5(l + ||c.||i.nc)^(A*:^)) UQ:P')\\ 
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Now by (|2.3ip applied with cr = 0, we have 



\D'^\^AM dry + aVl^)'ll^lli^nc^(M, 

Hence, for all n ^ 1 and for all x e {0, 1}"+^ such that a;o = 1, we have 
(2.43) 

o / \ n-\-l 

|||i?°|*G„,.||e. < (c(M)(l + «Vl^)'(l + ll^lli^ncMM,^)) m^p'Wx- 



It remains to treat the case a; 7^ (0, . . . , 0) with xq 7^ 1. Let n ^ 1, a; e {0, 1}"+^ 
such that a; 7^ (0, . . . , 0) and a;o = 0. Then there exists a smallest index jq > 1 such 
that Xj„ — 1. We can thus write 



« n 



3=ja + l 



Using (I2.32P jo times with = we get 



\d''\Ig. 



^ < (c(M)(l + aVT^)'(l + ll^lli2nc)) ''||(g,p')ll'^'x 

P n 

X J \\\D'\^A,{rj) n (^A.,W)l|e.d77. 

We estimate the right side using the estimate p.43p in the case xq — 1, and we 
obtain 

©2 ^ ^ 

/ \ n+1— Jo 

X (c(/i)(l + av/biAf(l + ||c.||i.ncWA^,^)) ||(Q,p')lir'"- 

Summing up, we have for all n 1 and all x, 
(2.44) 



|i?°|^G„,. < (G(A.)(l+aVl^)'(l + ll^lli^nc)^^(M,^))"^'ll(Q,p')llA 



We now estimate 



62 



and 



for all n > 1 and 



62 



X (0, . . . , 0). We treat both cases at the same time by estimating 

|||Z?TG„,.|Z?°r||e,, 

for any s,t E [0, 1]. For x 7^ (0, . . . , 0), there exists a smallest index ^ jo ^ n 
such that Xjg = 1. We distinguish three cases: jo = 0, jo = n, and 1 < jo < rt — 1. 
First, if jo = 0, then 

n-l 



|||i?"rG„,.|D"|*||g^^ / \\\DrAMX\\e.f[\\A.Av)X\ 
Writing 



(2.45) \D'rA,{ii)X 



D\ - pL + IT] 



ZJref (w) - ^i + irj 



-\n'\\D''\-'x 
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We obtain from the estimates above that 



\\{Q,p')\\x, 



and one concludes that 

/ \ n+l 

Then, if jq = one has 

|||7?°rG„.,|Z?°r|| / |||i?"|Mo(77)|U2^LHI^^o(^)|l24L^||^Ai(ry)|i^«nie.d77 

As before we have 



so that 



FinaUy, if 1 ^ jo =^ — 1 , then we write 

|||Z?°rG„,.|i?°r|L^ ^ / |l|i?"rAo(,7)llL^^LHl^^oWlliV4i.|lXAi(,y)X||e.x 



n-l 

X JJ P,^.(r?)X|U.^i2||A,„(77)|i?°|*|U2^i2d77, 
i=jo+i 



and using that 

\\XA,{ij)X\\e, 



c{^i){l + a,/W^f\\Lo\\l,^ce{^l,Lo) 



obtained by the same decomposition as in ()2.45p . we find 



/ \ n+l 

\\m'GnMD'\%^ ^ (G(M)(l + aybiA)'(l + ll^lli^ncMM,^)) \\{Q,p')\\ 
Hence, for all s, ^ G [0, 1], for all x 7^ (0, ... , 0) and all n ^ 1 we have 

\ 7l-\-l 

g(m)(i + aV^)Hi + Mhncn^^^^)) \\iQ,P)\ 



'Ml" 



'Ml" 



It remains to estimate 



G„,,|i?0|i 



The strategy is the same as the one for the 



estimation of 



, except that we first consider the case — 1. For 



any x G {0, 1}"+^ such that a;„ — 1, we have 
(2.46) G„,,|Z?"|i ^ f \f[\\A,,{v)X\\ 
In a similar fashion as the inequality (|2.42l) . we have 



A,ir^)\D^\l 



drj. 



l[\\A,^ifj)Xh.^L. ^ {Cif,)il + aVl^)Hi + MhncWiP^^y) \\iQ,P')\\x, 

3=0 
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62 ^ ( 







with now the factor d{fi,uj)'^ mstead of 9{fi,ijj) as m ()2.42p because of the term 
\\{Di-ci{uj) — + iv)^^\D'^\\\L2^L^- Using additionally (|2.33p . we get 

s ri+l 

+ a^l^y (1 + Mhnc) WiQ.p'Wx, 

for all X ^ (0, . . . , 0) such that Xn = 1- We treat the general case x ^ (0, . . . , 0) in 
the same fashion as the estimate of |i|-C)°|2Gji_2:||e2, using (I2.34p . and we get the 
same estimate as the x„ = 1 case. Finally, by (|2.29l) we conclude that for any n ^ 1 
and any x (0, . . . , 0), 

/ \ n+l 

(2.47) ||G„,,||gS^ (c(Ai)(l + aVl^)'(l + ll^llI^nc)^(M,c.)') WiQ^p'Wx- 

Combining this with p.27p . we can write 
(2.48) 

Vn ^ 2, \\GrAQ,p')\\Q ^ (c(m)(i + «v/b^)7 (1 + Ml^nc) SiP,^)')"^' \\iQ,p') 

while 
(2.49) 



\\GnQ,P)\\Q ^ [C{P){1 + av/logA)^ (1 + Mhnc) Oip,^)') v/logA||(Q,p')lk- 
We sum it up as 

(2.50) m{Q,p')\\ ^ «v/biAS||(Q,p')IU + E («2||(Q,p')ll 



n^2 

with 

E C(m) ((1 + «^/biA)' (1 + Mhnc) OiP,^?)^ ■ 

This is of course a rough estimate, but it has a simple form. We now turn to the 
estimates for the density terms, that is $2 • We start by estimating = po^ for 
n > 2. As for G^, we write 



ZTT ^ — ' 



2;e{0,l}" + ^ 

with Jn X '■= Pg„ X- It was proved in |10[ Proposition 10] that 
(2.51) 

Vn ^ 2, II J„,(o,....o)llL^nc {c{p){l + a^l^fa + ll^llW)) ll(Q,p')ll^, 

To estimate || ^n.^lU^nc for x (0, . . . , 0), we use the following lemma. 

Lemma 19. For any operator Q on i^A such that \D'^\^^°'Q G ©2 for some a > 
and such that l-D^^IQ G ©e/S; then pg G L'^nC and there exists a universal constant 
C such that 

(2.52) IIPQlU^nc ^ C (|||i?°|^+^Q||e. + |||i?"|g||6e 

Proof. We prove it by duality. Let V G Sja be a smooth function. Then, for any Q 
as in the statement of the lemma, we have QV G 61 with two different estimates 
on HQV^Ilei- First, we have 

IIQ^IIe. ^ ll|i?°|^+^QI|eJ||i?°r^-^^||e. ^ q||i?"|^+-Q||eJ|V^||L^, 
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by the Kato-Seiler-Simon inequality. This proves that pq £ and that 

IIpqIU. «;C|||Z?°|i+'^Q||e.. 

Secondly, we have 

\\QV\\e,^\\\D'\Q\\e.,MD^\-'V\\e, 

^ C\\\D^\Q\\e,JV\\L. < C\\mQ\\e,JS/VU., 

again by the Kato-Seiler-Simon inequality and the Sobolev inequality. This proves 
that PQ G C and that 

\\pq\\c^C\\\D^\ 



□ 



Let us first treat the case Xq = 1. Let n ^ 1 and x ^ (0, . . . , 0) with xq = 1. 
From inequalities (|2.4ip . (|2.3ip with any < cr < i, and (|2.35p . we find 



/ \ n+l 

\D'\^+^Gn.. ^(c{p)il + ay/]^)Hl + Ml..^cnp,io)) \\{Q,p')\ 



To deduce the general x case from the xo — 1 case we use as usual (|2.32p and (|2.36p . 
We thus have for all n ^ 2, 



n+l 



It remains to estimate \\Ji{Q,p')\\L^nc- We have proved that for x ^ (0,0), 

II Ji,.(Q,p')IU^nc < C(Ai) ((1 -t- a^f^)\\ -|- ||a;||i.nc)^(M, c^)) ' ll(Q,p')ll 
Let us write 



D\-p + tV 



so that 



-Rc 



-4oi(ry), 



27r ./» DI-P + i-q"^ D^- p + irj 



27r 7r J5a - M + 'i'V 



XAoiiri) dry 



+ P 



1 

2^ 



/ Aoi{rj)X ^ ^ . d7? 



Aoi(r/)XAoi(r/) d?7 



a:#(0,0) 



By [TOl Proposition 10], we deduce that 



IRr(Q,p')llL^nc ^ (C(Ai)(l + ayb^)^(l + ||c.||i.nc)e(/^,^)j ^J^m.p')\\x■ 
This ends the proof of Proposition [T71 □ 
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Proposition 20. Let v e L^(R^) n C and fi e (-1, 1). Let e > such that for all 
s G [1 — £, 1 + e], /i ^ as){D^'^). Assume furthermore that 

(2.53) / \og{2 + \k\f \V{k)\^ dk< +00. 

There exists ai = ai{^,v,e) > 0, Ai = Ki{^,v,e) ^ 1, and Li = Li(^,v,e) > 
such that for any ^ a ^ ai, for any A ^ Ai satisfying a log A ^ Li, and for any 

KG [1 — e, 1 + e] we have 

(2.54) II (Qrcf, Prcf)|U s^C(Ai, :/,£). 
Proof. Wc split Qrct into two terms, Qrcf = I + II, where 

I = X(-co,H(^r"')-^A.-+Gl^o(^^^re„), 
II = X(-oo,H (-Drcf) - X(-oo,H (-Da'''""). 

Since S L^(R'^) nC, we have Viy G and the eigenvalues of D^"^ are continuous 
with respect to s. Thus, there exists e' = s'{fi, i^,e) > e such that fi ^ (Tf,(D'"^) for 
all s e [1 — e', 1 + e']. Hence, there exists — Loil^TV^e) > such that for all 
KG [1 — £, 1 + e], for all a ^ and for all A ^ 1 satisfying a log A ^ Lq we have 



K 



e [l-e',l + e']. 



1 + aBA 

Using Lemma fT2l this implies d{fi, a^,^ (/jKi/^n-j-j ^ <^o(/^, £, io(M, i', e)) > 0, for all 
KG [1 — e, 1 + e], for all A ^ Ao(/x, z^, e, Lo(Mj i': e)), and for all a ^ such that 
a log A ^ io- In particular, ^ asjj^{D'^'^"'") and we can write 



27r DO - ^ + 177 £'a - + ^r'°" -t^ + ir]' 

As in the proof of Proposition 1171 we estimate this term by decomposing it into 
L = Ii + I2, with 



^1 = ^ / f no \ ■ nAK.,.„nA) X 



2 

X 

J=0 



^2 = Tr / ( no \ ■ nAK...„nA ) ^ ^ . dr/. 

By [101 Lemma 13, Lemma 15], we have 

||/l||QS^C(l + ||«J.,e„||i.nc)- 

Notice that 

27r Jr 

where ^1(77) was defined in (H^H). By (H^H), (lOTI) . and the trivial bound 



2|||i^Vi('?)|a=o,™„Ji^"l'/'l|e.+2|||i^T/'^i('7)|a=o,.=..„Ji?"| 



©2 



_^ C(M)lk^'rcn||i2nC^(A^^ ^^'rciQla^O 
£^(77)3/2 
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we infer that 

II^2||q C{^I)\\Kl',cn\\hnC^i^^,^')fa=0■ 

Hence, for any k e [1 — e, 1 + e], 

||/||q ^ C(/.)(l + (1 + £)'|l^^|lI.nc)^(M,«^rc„)f„=o- 

Since d(/j,,(Tj5^ ^ ^o, we have ^(/i, Ki'i.cn)|Q=o ^ C{^,v,e), uniformly in 

a ^ (the dependence on a only being in Vrcn, which is not obvious from the 
notation), A > Aq satisfying a log A ^ Lq. This leads to 

\\I\\Q^^C{^i,v,e). 

By (|2.10l) , this implies in particular that for A ^ 3 so that log A ^ 1 , 

(2.55) \\X(-ooM {DT-)-PI_\\q ^ C(M,;^,e)v/b^. 

Using again Lemma [121 we know that for all ^ a ^ ao and for all A ^ Aq 
satisfying alogA ^ Lqj we have d{fj.,asjj^{Drcf)) > Co for aU k G [1 — £, 1 + e]. 
Hence, we can expand II by Kato's formula, as for /. To estimate this expansion 
by the same method as in the proof of Proposition [TTl we need the additional 
condition 

8a||Kfrcn||c < 1 - ImIi 

which ensures the finitess of 9{fi, nv-^cn), as seen by Lemma 1131 This condition is 
satisfied for all ^ a ^ ai, for all A ^ Aq satisfying a log A ^ Lq, and for all 
K g [1 — e, 1 + e], for a certain ^ ai = ai(/i, v, e) ^ ao(^, v, s). Using the same 
method used to estimate ||/||q and the proof of Proposition [T7l we find 

(2.56) ||//||q 5$ C(M)a\/k^(l + a^k^^ {1 + (1 + ef Ml2nc)^{ti, ni^rcnf . 

We still have 9{iJL,Kv^cn) ^ C{ii,v,e) for all k e [1 — £,1 + e], and we also have 
aVlog A ^ a log A ^ io if A ^ 3. Hence, 

\\II\\Q^C{ii,y,e) 

for all ^ a ^ ai, for all A ^ max(3,Ao) with a log A ^ Lq and for all k £ 
[1 — £, 1 + e]. The estimate on prcf is more subtle. We decompose it as 

Prof -UAKUicn + p [I] + p [II] ■ 

The operator Ua is the multiplication operator in Fourier space by the function 
[/A(|fc|) := Ba — BA{k). It has been defined and studied in [7]. In fact, combining 
[71 Equation (3.3)] and 7, Equation (A. 3)], one gets the estimate 

(2.57) VA^l,Vr^O, ^ UA{r) ^ m (^l + ^\og{l + r^)^ , 
with Ki — 258/7r. We deduce 

II^AKi^rcnllc =^ C(l + e) |1 :/|| i2nc , 

\\UAniyren\\l2 ^C{l + ef f log(2+|fc|)2|P(fc)|2dfc. 
Next, using the same ideas as before, one gets 

IIpWIIl^CC C{p){l + (1 + sfMlh^cWip, KU,en)l=o ^ C{p, V, e), 

\\p[II]\\L-nc ^ C(M)alogA(l+a\/biA)^(l+(l+£)^||i.||[.nc)^(A'.'^'^re„)' < C{p,i^,e), 
which ends the proof. □ 
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2.1.4. End of the proof of Theorem 1. We apply the fixed point theorem of Banach- 
Picard to the function $ defined earher by 

for any {Q,p') G Bx{R) C X, where 

Bx{R) = {(g, p) G X, II (Q + Q,et, < R], 

for a, R > to be determined later on. For ^{Q,p') to be correctly defined, 
we have to check the assumptions of Proposition [iTl First, as in the proof of 
Proposition [20l using Lemma [T2| there exists Lq = Lq{ijl,v,£) > such that for 
all ^ a ^ ao{fi,i',e, Lq), for all A ^ Ao(/i, z^, £, £o) with alogA ^ Lq, we have 
p ^ iTj--i^(I?iof) for all K G [1 — £, 1 + e]. Next, we assume additionally that 

8a||Kz/rcn||c < 1 - \p\, ab\\{Q + Qrc{,p')\\x < 1- 

This condition is verified if we assume furthermore that ^ a ^ a2 for a certain 
Q!2 = Oi2{p,v,e,R) > 0. Using Lemma[T4l it implies that /i ^ <Js^f^{DQ), where we 
recall that 

= X(-oo.m] (Acf + anA(Vp' - i?Q+Q,,^j)nA). 
We can now use Proposition [TTl to infer that for all (Q, p') G Bx{R), 

(2.58) ||$'(Q,p')IU^<=^ ^«\/l^S + a^n(aSi?)"-i. 

Assuming A ^ 3, we have a^log A ^ a log A ^ Lq. Moreover, as in the proof of 
Proposition [20l we can show, using Lemma[T2l that there exists C — C{p,v,£) > 
such that for all ^ a ^ uq, for all A ^ Aq with a log A ^ Lq, and for all 
KG [1 — e, 1 + e] we have S ^ C. As a consequence, the right side of (|2.58l) is finite 
as soon as a ^ {CR)~^ /2 = = a^lfj,, v, e, R). It thus satisfies 

V(Q,p') G Bx{R), W{Q,p')\\x ^ C,^, 

for some Ci = Ci{p,v,e,R) > 0. Hence, $ is a contraction on Bx{R) if a ^ 
(4Ci)^^ = = ai(p,v,e,R) > 0. It remains to prove that Bx{R) is invariant 
under $. For any {Q,p') G Bx{R), we have 

\\HQ,p) + (Qrct,o)||;, WHQ^p) ~ H-Q..iM\x + ||$(-Qrcf,o) + (grcf,o)|U 

^ Cl + I KQrof , Prcf ) 1 1 A' • 

By Proposition [201 if ^ a ^ ai, A ^ Ai, and a log A ^ £i, we have 

||(Qrcf,/5rcf)|U ^ C ^ C {p, V , s) . 

Since we have chosen a such that Ci^/a ^ 1/2, $ stabilizes Bx{R) if 

R = 2C. 

To recap, if ^ a ^ min(ao, ai, a2, as, 04), A ^ max(3,Ao,Ai) with alogA ^ 
min(io,-^i), then for any k G [1 — £, 1 + e], $ is a contraction on Bx{R), with 
R = R{p, V, e) > 0. Hence, there exists a unique Q £ X satisfying 

Q X{-oo,p.] (^rcf + allA [VpQ+p^^f ~ i?Q+Q„f] Ha) - Pref, 
||(g + Qrcf J PQ Prcf 
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Defining Q := Q + Qrcf - G'i,o(Ki'rcn), we iiave pg + prcf ^ Pq ~ B^Kv^cn and we 
tiave proved that tliere is a unique solution to 



Q^X{-oo,tA Lor + "Ha 



Ha - Pa" 



A,- 



II (Q + Glfi{KV^cn), Pq - BAKV^on)\\x ^ R, 

which ends the proof of Theorem [TJ 

2.2. Proof of Theorem m Let L2(R3)nC, p± G (-1,1), and £ > satisfying 
the assumptions given in Section 11.21 Since the eigenvalues of D^^ are continuous 
with respect to s, there exists e' > e such that (/i±, v, e') also satisfy the assumptions 
given in Section [TT^ By the same method as in the proofs of Lemma [T^ and Propo- 
sition[20l there exists a' — a'{p,i',e) > 0, A' = K'{p,,v,e) > 1, L' = L'{p,,i/,e) > 
and C — ^'(/i, I', e) > such that for all ^ a ^ a', for all A ^ A' satisfying 
alogA < L' , and for aU k e [1 — e, 1 + e], we have k/(1 + aB^) e [1 — e', 1 + e'] 
and 17^5(1?^'^"")) ^ The following result ensures that the operator D'^'°'' 

also has an eigenvalue in for A large enough. 

Lemma 21. Let v ^ L"^ C\C, p± G (—1,1), and e > satisfying the assump- 
tions given in Section \l.2l Then, there exists A" = A"{fi±,i>,e) ^ A' and C = 
C{p.±,i^,s) > such that for any A ^ A", for any ^ a ^ a' with alogA ^ L' , 
and for any k G [1— £, 1+e], the operator D']:^""^ onijA possesses a simple eigenvalue 
A(Ki^rcn,A) in {p^,fi-f-), which verifies the estimate 

(2.59) |A(Ki.,en) - A(«;i/„„, A)| «C CE{A)-\ 

Furthermore, if we denote by (/^(Kt'ren, A) the unique normalized eigenvector of 
-Da""^™ for the eigenvalue A(«;fron,A) such that ((p(Ki/rcn), 95('«^'ron, A))i2 > 0, then 
we have 

(2.60) MnVren) " V?(Kl'rcn , A) 1 1 < CE{A)-\ 

Remark 22. Since k/ (1 + uBa) G [1 — e', 1 + e'], and {p±,i', e') satisfy the assump- 
tions of Section [L2l A(Kt'ron) = -^(«^/(l + ctBA)i^) is well-defined. 
Proof. In order to compare D'"^"" with D'^"'"" , we identify 0'^'°" with the operator 
on which is D'^""" on S^a and on (1 — nA)iOA- Let "^(o be a complex, compact 
contour around {p.^, P+), of index 1, crossing R only through /i_ and This way 
'^^ remains uniformly far from ctj, (i?'"'"" ) and (Tj-,(£'^'^'°"), for A ^ A'. We then 
have 

A(M-.M+)(^'''^"") - nAX(M-,M+)(^r™)nA = 



^ ^nA(— iHAdz 



2iTT Jr^ yo'^"'-" - z D'J^"'' - z 
" i L " ^^^I^^"^ + nA^^;^(l - Ha) 

+(1 - nA) j;..,i,, a - Ha)) dz ^ / + //, 
where / is only the term in the second line. For the second term, we write 



\II\\l-^L-2 ^ 3||(l-nA)|i?'riL^^L^SUp 

C{p±,v,e)E{K)-\ 
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The first term is computed by the resolvent formula 



- — -(1 - nA)K._nA-— i;- 



so that we have the estimation 



s$ sup 



|D"|-l(l-nA)|U2^i2 sup 



1 



As a consequence, for A large enough such that for instance 

-nAX(,,_,^+)(i?r™)nAlU^^L^ < 1, 

the operator D'^""'" on Sja has a simple eigenvalue A(Ki'rcn,A) e To 
estimate the difference between these eigenvalues, notice that we have proved that 
for any z G [/i± — <5o, M± + '^o], for some = So{^± , i', e) > small enough 



■f3A-i--f>A 



1 



-Ha 



^:C{^i±,v,e)E{K)-\ 



It implies that, for A large enough, 

1 1 



X{KUren) - Z X{KVrenj A) - Z 

and hence 

|A(Kfrc„) - A(Kf,.o„, A)| s; C{ii±,v,e)E{K)-'^. 
Let us now consider the unique eigenvector ip{Hiv-ccm A) in the statement of the 
lemma. It is well-defined for A large enough since the operator 

A := \Lp{KV^en)){'P{KVren)\ " | (/^(Kl'ren, A)) ((/3(Kt'ren , A) | . 

verifies 

^ {Aip{KV^^n),Lp{K.V^a^))L2 = 1 - ((y5(Kl'ren, A),93(Ki^rcn))i2 

^ \\A\\l2^l- ^ C{fi±,iy,e)E{A)-' < 1. 

Hence, {(p{KVrcn, ■^), v{nvicn)) 7^ and 

||</?(Kt'ron) - V3(Kt'ren, A)||^2 = 2(1 - {ip^KVrcn) , ^{Kl^Tcn, A))) 

SC 2(1- ((p(KZ.,en),(^(Kl'rcn,A))2) <; C {^l± , V, s) E{Ay\ 

To obtain a control in H^, we use the eigenvalue equation, which implies 

D°{(p{KVrcn) - (p{KVrcn,A)) = A(KZ^i.on) (<p(Ki^rcn) - (y3(KJ^rcn, A)) 
+ (A(Kt'ron) - A(K^',.on, A))(^(Kt'rcn, A) 
+ VKi.„„((y5(Ki^ren) - <y3(Kl'rcn, A)) + (1 - IIa) Kz.^,, (/^(Kl'i.cn , A) . 

Using the fact that VKi,^^^ip{Kh'rcn, A) is bounded in i/^ uniformly in A, we get the 
estimate in iJ^. □ 

We now begin the proof of Theorem [21 Since we work with k fixed we will 
omit the dependence on k in our notations, meaning that we write Q±{a) instead 
of Q±{K,a). We assume ^ a ^ min(ao,a')j ^ ^ max(Ao,A") and alogA < 
min(Lo,L'), where (ao,Ao,-Lo) is given by Theorem [T] and (a',A",L') were given 
just above. In the sequel, we will use notations such as Oa^o{c() (resp. O^^Q{a)) 
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to denote a function of a which is bounded in absolute value (resp. in norm || • ||x 
for some space X) by Ca, where C only depends on /i±, e. In our regime of pa- 
rameters, notice that we have for instance Oa-^ooCVIoS^) = Oa^Q{\foL). Likewise, 
in the following we will use the notation C for a constant only depending on [i± , v 
and £, whose value may change from line to line. 
We use the a priori bound given by Theorem [T] 

(2.61) ||(Q±(a) + Gl,o(Kl'rcn),P±(a) - -BAKt'rcn)|U = Oa^oll), 

where we used the shortcut notation p±(oi) := pq^^a)- The self consistent equation 
for Q±{a) is 

Q±ia) - X(-oo.^±] (nA(i?r + «(V;±(a)-i?Q±(a)))nA) -P?,- 

= X(-00,AI±] (A-Cf + (Vp^(Q)-SAKl'ro„ ^ RQ±(a)+Gl,o{K,^,an)) ^^a) - ^A, 

= ^ a"'Gn,± {Q±{a) + Gifi{KI^ren),P±{a) - BAKI^ren) + Qref,±, 

where we used a notation similar to Section 12.11 (with a slight modification of the 
definition of Qrcf,±)- 

Qrcf.± ■= X(-oo,M±] (A-cf) - Pa-, 

D,,f nA(i?"''"" + ai?i,o(KJ^ro„))nA, 

Gn,±iQ,P)--=^-^ / 7^ --[(^p-^Q^n T-) 

27r Jm Dy-ct - p,± + 17] \ D^cf-p±+ivJ 

Notice that Gi.o(Ki'rcn) is independent of the index ±. By (|2.48p . we have for all 

\\Gn,±{Q,P)\\Q^C\\{Q,p)r, 

so that 

\\Gn,±iQ±{a) + Gifl{KVrcn),P±ia) ~ BAKiy,cn)\\Q ^ C" . 

On the other hand, 



with 



Gi,±(Q,p) = Gi,o^±(p)-Go,i^±(Q), 

Gi,o,±(p) ■=1^ [ -R ^——TIaVpUa- ^— -dr?, 

27r iJi-cf ~ P±+iri Diof - M± + 

Go,i.±(Q) -^l^ [ 75 ^-^UaRqU^- ^-^dr]. 

27r Drcf - fJ.±+ nj Act - ^^±+^1^ 

Combining p.47p and [TO^ Lemma 11], we get 

l|Go,i,±(g)||Q ^ GIIQIIq, ||Gi,o,±(p)||q ^ GV^iAllplU^nc- 
Hence, using (|2.6ip we get 

\\Q±{a) - Qrcf,±||g = Oa^a{^/a). 

Now using (|2.56p . we deduce 

(2.62) |i(3±(a) — Qlin,±(KJ^rcn)||g ~ Ca-i-0 

where 

Qlm,±(^)=X(-oo,M±] Pa)-^A.-, 
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which in particular implies Corollary |3l The self-consistent equation satisfied by 

p±{a) is 

P±(a) = Prcf,± - aBA{p±{a) ~ BAnv^cn) + aUA{p±{a) - BAKv^cn) 

+ aPl,0,±,r(p±(a) - BAUVrcn) ~ aPoa,± (Qi («) + Gifl{KU^cn)) 

+ ^ Pn,±{Q±{a) + Gi,o(Kt'ron),P±(a) - BAKVrcn), 

where p„,±(Q,p) := Pg„.±(q,p) for all n > 2, pi.cf,± := PQ„f,± and for any p, Q, 



Pi,o,±Ap) P 



1 



1 



2n 7r Dy-cf - p±+iri 



UaVpUa 



1 



1 



-Drcf - P±+iv 



dry 



^ nAT^pHA ^„ d77 



Poa,±(Q) = P 



1 



27r 7k + iV D° + irj 

1 



—UaRqUa- 



_27r 7m Acf- M± " Dj-ef - p± + iri 

We rewrite the self-consistent equation for p± (a) as 

(2.63) 

P±(a) - BAKV^cn = -r- 5-(/0rof,± - BaKV^cu) + — 5-Wa(/0±(q:) - BAKV-ccn) 



I + uBa 



1 -I- a^A 



l + a^A 1 + q;Ba 

+ , ^ „ Q:"yOn,±(<3±(a) + Gi,o(Kfc'ron),/3±(a) " BAKV^cn)- 

1 + aSA ^ 



n^2 

We now have the necessary estimates to compute the limit of F{k, a) as a — > 0. 
2.2.1. Kinetic energy terms. By (j2.62p . we have 

(2.64) Q+{a) = Q-{a) + \Anu,,^, K)) {^{kv,,,,, K)\ + O^^^i^) . 
For any projector P such that Q = P — P'^ E ICa we have 

(2.65) Tro(Z?"Q) = |||i?°|^Q|||, ^IIQIII. 

The subtlety of this estimate is that it only controls the kinetic energy by the Q- 
norm only for Q's that are of the form P — P^. Indeed, it would be wrong to deduce 
from that 

Tro(i?"Q+(a)) -Tro(i?"Q_(a)) = {^{nu,,^, K) , D\{nu,,^, K)) + O 

since the operator Q^{a) — Q-{a) — \ip{KVT:c-a, f^)) {fii^v-ccn, ^)\ is not of the form 
P — po fQj. some projector P. Instead, one should see Tro(I?°Q±(a)) as a quadratic 
term by the first equality in (|2.65l) . Writing Q±{a) = Qim,±(«^^'rcn) + Q±,r{ct), we 
thus have 

(2.66) Tr(|i?"|Q±(a)2) = Tr(|i?"|Qn„,±(«i/,e„)') + Tr(|Z?"|Q±,,(a)2) 

+ Tr(|i?0|{Qii„,±(Ki/rc„),Q±,r(«)}), 
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where {A, B} := AB + BA. The first term has a good behaviour since (3iin,±('«^'rcn) 
is of the form P — P*^, hence the first term is linear in (3iin,±(Ki'icn) and we have 

(2.67) Tr(|Z?"|Qun,+ («^^rcn)') - Tr(|i5"|Qun,-(«^^rcn)') 

= TTo{D"{Qiin,+ iKiyrcn) - Qlin - ('«^'ron))) 

= (i^(Kt'rcn, A),i:)°<^(Kt'rcn, A)). 

The second term in (12.66^ is controlled using that ||(5±,r(a)||Q — Oa^o{y/a), hence 
Tr(|i:'°|Q±.r(a)^) s$ ||(3±,r(a)|| q = Oa^oCa)- Finally, the term with the Poisson 
bracket in (|2.66p cannot be shown to go to zero as a — >■ 0. Indeed, it is just a 
Oq^o(1) since ||(3iin,±(Ki^rcn)||Q = CA->oo(VlogA) by p.55p and ||(5±,r(a)||Q = 
Oa^o{^/a). We will compute explicitly the non- vanishing term in this expression. 
We have 

(2.68) Tro(i?°Q+(a)) - Tro(i?"Q-(a)) = ( ^{nUrcn, A), D'^ip{Kiy,.cn, A)) 

+ Tr(|Z?°|{Qn„,+(Kz.,e„), (Q+,.(a) - Q-Aa))}) 

A), Q_,,(a)}¥'(«:i/rcn, A)) + Oa^o{a). 

The computations are lengthy and done in Appendix [A] We obtain: 

(2.69) Tro(i?"g+(a)) - Tro(i?"g- (a)) = ((^(Ki^.-en, A), i?V('«'^rc„, A)) 

^ Q^^A ^y)(Kt^rcn, A), VKi/„„y(Ati/i.cn, A)) +e'a^o(a)- 

1 + aBA 

In particular, notice that 

lim Tro(I?°Q+(a)) - Tro(I?°Q-(a)) ^ {ip{Kiy,,^) , D\{Kiy,,^)) , 

a— >-0 

meaning that the difference of kinetic energy of the interacting vacuums does not 
converge to the difference of kinetic energy of the non-interacting vacuum. This 
shows the subtlety of the limit a — >■ 0. Fortunately, the extra term in ()2.69p will be 
compensated by another one coming from the direct term. 

2.2.2. Potential energy term. The second term to compute in F is the only "true" 
linear term: that is ~D(ki^, •). By (lA.ip . we have 

(2.70) p+(«) -p_(a) = J |y.(^t.,en,A)p + 0^!ff(a), 

1 + aBA 

which implies that 

~ D{Ki^,p+{a) - p^{a)) = -— — -L»(ki^, |93(Kz^ren, A)p) + C'„^o(a) 

1 + uBa 

(2.71) = -((^(KZ^rcn, A), VKi,^^„(p(Kt'rcn, A)) + C'Q^o(a)- 

2.2.3. Direct term. The next term we have to treat is 

^(^D{p+{a))-D{pAa))) = |(i^(p+(a) - p_ (a)) + 2I?(p+(a) - p_ (a), /,_(«))) , 

where D{p) is a shortcut notation for D{p,p). By (jA.ip and (|2.6ip . we have 

aD{p+{a) - p^{a)) = Oa^a{a), 

aD{p+{a) - p^{a),p^{a)) = {A^v^cn, J^) iV^u„,„v{i^v^cni ^)) +Oa^a{a). 

1 + aBA 
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Hence, 
(2.72) 

2 V / 1 + aisA 

which exactly compensates the spurious term in (j2.69p . 

2.2.4. Exchange term. Finally, we treat the last term of the energy, the exchange 
term. For any operator Q, we have 

jR^xR^ \^-y\ Jr^xR^ \ E[p-\-q)2 J 

1 

< CWQWl, 

where in the last inequality we used [lO) Lemma 8]. The term we have to estimate 
is 



a 



dx dy 



2 VVr6 \x~y\ Jr6 \x-y\ 

a f |[Q+(a)-Q_(a)](x,y)|2 



2 Jr6 \x - y\ 



I Trc4 ([Q+(a) - Q-(a)](x, y)[Q.{a)]{x, y)*) 

\x-y\ 



By ([2:641) together with |lQ_(a)||g = O(VlogA), we have 
a f |[Q+(a)-g_(a)](x,2/)|2 



\x - y\ 



dxdy = 0^^0(0:), 



(2.73) Trc4 {[Q+{a) - Q-{a)]{x,y)[Q-{a)]{x,yr ) ^^^^ 

Jr^ \x - y\ 

= Q;((/3(KJ/ren, A),i?Q_(a).^(Kf,.on, A)) +Oa^Qia) 

Using that 



|((/3(Kl'ren, A),i?Q_(c,)V3(Kl't.en, A))| ^ |1 i?Q_ (q) || //i^L^ |i ¥'('«^'rcn, A) | 



2 



and the fact that \\RQ_{a)\\m^L^ = pQ_(a)l^"l ^IIl^^l^ ^ C||Q_(a)||62 
Cq^o(1): "we obtain 



Regrouping (P^Ml) . (PTT]) . ((^T^ . and together, we get 

F{k, a) = A(Kt'rcn, A) + Oa^oia). 
Since A(Ki/rc„) = A(Kl/i.on, A) + 0(£'(A)-1) = A(Kt'ren, A) + ©(e^i) = A(Kl'ron, A) + 

0{a), we deduce the theorem. □ 
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Appendix A. Kinetic energy estimates 
We expand Q+^r{a) ~ Q-,r{a) up to Oa^o{a'^^^)- Recall that we have 
Q±{a) = Qrcf,± + aGi,o,±(/5±(a) - Bahv^cu) 

~ «Go,i,±(Q±(a) + Gi,o(Ati^rc„)) + O^-^o(a')- 
Then, using similar methods as in the proof of Propositions [T7] and [501 we obtain 

Gifi,±{p±{a) - BAKt-rcn) = G^^d _t(p±(a) - Bf^nv^cn) + O^^Q{^/a), 

Go,i,±(Q±(a) +Gi,o(Ki^rcn)) = G[,';i'_±(g±(a) + Gi^o('«;^rc„)) + 0?_^o(V^): 
with the notation 

G[t±iP) -^ir [ ly^^ — nAV-pHA ^ — d,y, 

G^1±(Q) [ ^ — HaAqHa ^,,.,„„ ^ — dr;. 

27r - M± + ^a - M± + 

Hence, 

Q±(a) = Qrcf,± + aG^*}] j_(p±(a) - Bakv^c-o) 

- aG[%^{Q±{a) + Gi,o(Kz^rcn)) + O^_^o(«'^')- 

This leads to the decomposition 
with 

Ql Qrcf,+ — Qlin,+ (K^'rcn) — (Qrcf;- — Qlin,- ('«^'rcn)), 
Q2 ■= Gl^O,+ iP+i(^) ^ B^KV^cu) - Gi*}] - SAKt'ron), 

Q3 G<,;^i\+(Q+(a) + Gi,o(Ati/rc„)) - G^^l_{Q.{a) + Gi,o(Ki^,en)). 
We will treat each of these three terms separately. 

A.l. Estimate on Qi. To expand Qi, we write 

'5rof,±-(9lin,±('«t'rcn) = aG^j'^j' (Gi,o (KJ^ron) ) -a^G^j^^^ ^_ (Gi,o (Kl'rcn) ) +C'q_>o ("^^^ ) ) 

where for any Q we have 



2 

dyy. 



Hence, 

Qi = «(G^:i^-G^';i^,_)(Gi^o(A^^^rc„))"«2(G^«)^^^G(l_)(Gi,o(«J^rc„) 
For p = 1, 2 and Q G Q, we notice that 

(Git. ~ Gii,_)iQ) = ±- £ (nA^QHA^^^)'' dz, 

where is a contour around in the complex plane which intersects the real axis 
only at fi^ and By the residuum formula, one computes that 

(G^l+-G^l_)(Gi,o(/ti^rc„)) = |¥'( 

A))(5(K:^ren)v('«^'ron, A)| +C.C., 
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+ |S'(KI/rc„)<p(KJ^rcn, A))(S'(KZyren)</'(Kl'ren, A)|, 

where c.c. denotes the adjoint of the preceding operator and 

S{KV,cn) ■■= -f^;^^ ^, rTi?l,o(Kl'rcn), 

Df^ - A(KJ/ren, A) 

the operator P;^ being the projection on the orthogonal to <y9(Kfc'icn, A) in 55a- By 
(PT7)) . we have 

||5'(Kfi.cn)||//i^_H-i ==00^0(1)- 

Since the we have to estimate the kinetic energy of these quantities, we introduce 
a new norm corresponding exactly to the kinetic energy, for which the rank one 
terms above are easier to estimate. For any Q such that the following expression is 
well-defined, we set ||Q||kin := || |-D°|^/^Q||e2- We thus see that for any G 
lll''/')(Clllkin IIV'll_f/i||CllHi- Foi" this norm, we can write 

Ql = Q:(|99(Ki^i.on, A))(5(Kt'ron)'^('«^'rcn, A)| + 1 5(Kt'ron)'^('«^'ron , A)) ((^(Ki^rcn , A) | ) 

Using now that Qiin.+ ('tt'rcn) = ^ Gi^o('*^'rcn) + Cq-^o(-'^)' deduce 

Tr(|7?°|{Qun,+ ('ti^rcn),Ql}) 

= -2Q!Re(i^(/tt'rcn, A), {\D'^\,Gifl{KV^cn)}S{KV^cn)'P{nV^cn,^)) + O a-^Q{a) . 

Next, we have 

I ((/^(Kl'ren, A), G'i,o(Kt'ron)}'5'(Kt'rcn)</?(KJ^rcn, A)) I 

^ 2||Gi,o(Kfc'rcn)||e2ll'S'(Kl'rcn)|l_H-i-+L2|l</'('^'^ron, A)||^i = O a^Q{l) , 

SO that 

lin,+ (/^^rcn ),Qi}) =Oa^o(a). 
A. 2. Estimate on Q2- We split Q2 as 

In the same fashion as before, the residuum formula leads to 

|(p(Kt'rcn, A))(S'(V+).^(Kt'rcn, A)| + 1 )<^(Ki/rc„ , A)) ((^(KZ^ren , A) | , 

with V+ := V^p+(q)-Baki/,,„ and 

Notice that we have ||S'(F+)||i2_j.^i = Oa^Q{l). The second term in the expansion 
of Q2 requires to develop p+{a) — p-{a) up to 0{a). Combining p.63p with the 
estimates valid for any p € n C, any Q G Q, and any n 5^ 2, 

\\Pn,±{Q,p)\\L^nc < G\\{Q,p)\\, \\pifl,±Ap)\\L^nc < ^11^11^2^,^, 

P0,1,±{Q±{ol) + Gifi{KV^cn)) = P0,l((3±(a) + Gifl{KV^ 

en 
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where po,iiQ) '■— p[G'o,i(*3)]j deduce that 
1 



p+(a) - p-{a) 



a 



1 , n (Prcf,+ - Prof,-) + 

1 + aB\ 1 + aBA 



Po.i{Q+{a)-Q-{a))+O^^^S{a) 



1 + aBiy 



Since /9i.cf,+ - Prcf- = |(/3(Ki/rcn, A)p + pQ^ and 



Pq, = 2aRe (^.^(Kt'ron, A)S'(KZ^,.cn)v('«^'rcn, A) j + O^^O 

it remains to estimate ||Re(V'C)l|L2nc ^oi any ■0,C G First, by the Sobolev 

imbedding ^ we have |lRe(V^C)l|L= «J IIV-IIl^IICIIl^ < Mh4C\\h^- Sec- 
ondly, by the Hardy-Littlewood-Sobolev inequahty, we have 



Hence, 



Re (/'('t^'ron, A)5'(Kfren)<P(Kfren, A) 



Oa^o(l), 



and thus Prcf,+ ^ Prof,- = |v'('*^ron, A)p + O^^^f («). For the third term in the 
expansion of p+{a) — p-{a), we use Q+(a) — Q-{a) — \^p){(p\ + 0^{y/a) and the 
estimate \\pa,iiQ)\\L^nc ^ C^log A||Q|| g to infer 

po,i(Q+(") - - poA\'p){v>\) + o^^Sii). 

Lemma 23. We have the estimate 

IIPo.id'/'X'y'DllL^nc = Ca^o(l)- 

Proof of LemmaW^ We estimate the Coulomb norm by duality: for C G C' a 
smooth function, we have 

|Tr(Go,i(l¥')('^l)C)l - |Tr(i?|^)(^|Gi.o(C))| P|^>(^|ll6j|Gi,o(C)l|6. C||<^||^, ||C||c'. 

Hence, |ipo,i(IV')('^l)llc — Oa^o[^). Wc cannot use the same method to estimate 
the L^-norm, for which we have to be more precise. Let us use momentarily the 
shortcut notation po.i Po,i(|'/')(vl) f-i^d recall that for any fc S we have 



|^+fe/2|sJA 



Trc4 - fc/2, 1 + fc/2)M(£ - fc/2, i + fc/2) 



^|^>(^|(^-fc/2,^ + fc/2) 



< C\kY i?|^)(^|(^-fc/2,£ + fc/2) 



|Af(^-fc/2,£+fc/2)|2 



d^, 



where we used the estimate |M(p, g)p ^ C\p — q\'^E{p + (7) * (recall that M(p, 5) 
was defined in (|2.13p '). As a consequence, we have 



Wpq.iWl^ ^ C \p - q\ dpdq = C|l[-iV,i?|^)(^|]|l(e^)3. 



But since [-iV, ?/) 



Vy(a:)y(y)* + y(3:)Vy(j/)* 



we infer from the Hardy 



inequality that Upo.iIIl^ < C\\ip\\jji, so that Upo.ilU^ = Oq^o(1)- 



□ 
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The last term to estimate in p+ (a) — p_ (a) is Z^a (p+ (a) — p- (a) ) . Using (|2.63p , 
we have for all fc e M"^, 

UA(\k\) 
1 + aBA(k) 

' "^aC'^I) ■F[pi,o,±,,.(p±(a)-i?AA.t/,en)](fc)- "^^yP,, J-[po,l,±(Q±(a)+Gl.o(^i^rcn))](fe) 



1 + aBAik) ^-^."'-''^--^"^ 1 + aBAik) 

, aUA{\k\) 



1 , D ^J'[p„,±(Q±(a) + G'i,o(Kl'rcn),P±(a) - BAKZyrcn)](fc)- 

1 + aBA[k) ^ 



Since 



aUA[\k\) 



l + aBA{k) 



aB 



A 



aBA^k) 



1 + aBAik) 



we have 

UA{p±{a)--BAKV^cu) = VA(prcf,±-5AKi^rcn)-aVAPoa(Q±(a)+Gi,o(Kfc'rcn))+C'^J^o'^(l), 

where Va is the Fourier multiplier by fc i— Y^^^^p;y- Hence, for the difference of 
densities we obtain 

Ua{p+{c^) - = VA(prcf,+ - Prct,-) - aVAPo.i (Q+ (") - + Oa^o'^(l)- 

We have already seen that po.ilQ+l^) ^ Q-(Q!)) = C'^!^g''(l), and that Pref,+ — 
Prcf - = [(^(K^'rcn, A)p + 0^4i'^(a), which implies 

UA{p+{a) - p-{a)) = VA|(^(Ki^rc„, A)|2 + 

By (P?57)l and the estimate (1 + aBA{k)y^ 1 following from BA{k) ^ 0, it is 
enough to know that |</jp is bounded in to prove that Z^a(p+(q;) ^ P-(q^)) — 
0^^q{\). Since ^ C||(/7||^2 and is an eigenvector of D'^^"" , it is suffi- 

cient to control || V(VKiyrcn'/')llL2- By the Hardy-Littlewood-Sobolev together with 
|VK| s$ Civ] -k I • v e implies that VK e i*^, and ip e ^ implies 
that (pV14,_ e Hence, (p = Of!,o(l) and UA{p+{a) - - 0^!ilf (1). 

After these painful estimates, we obtain 

(A.l) p+ia)~p^ia) = -_i_|^(^j.,,„,A)p 

which in turn implies 

A)| + |5(V+)(/3(Kl',.cn, A))((p(kI/,.o„, A)| 

This leads to the expansion 

Tr(|i?0|{Qii„,+ («:!.„„), aga}) = -2aRe((p(Kz/„„, A), Gi,o(KZ/re„)}5(l^+)(p(KZ/rcn, A)) 

+ " Tr(|i^°|{Qun,+ (^^rcn),Gl1,(|y'(A^^rcn,A)n}) +0„^O(«) 

1 + aisA ' ' 

As before we have 

\{^{>iI^ron,A),{\D°\,Gi,o{Kl^ron)}S{V+)ip{K,^ron,A)) \ = Oa^o(l)- 
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Furthermore, we write (3iin,+ (KJ^rcn) = -Gi,o(Ki^rcn)+Qiin,+,r(Kt'rcn) and G^'^o _ = 
G,M^) + g[^1^^,M^), so that 

+ Tr(|I?"|{gn„, + ,.(A£^^ren),G(;'),_,,(|^n}). 

We know that \\G^iI^^^{\^\'')\\q G|| |(p|2||^.nc ^ C\Mm ^ 0^^a{l), and 

||Qlin, + ,r(Kt'rcn)||Q = Oa^oCl). HenCe, 

Tr(|i?0|{Qn„,+,,(Kz.,en),Gl;=oV,,(|(pn}) = O„^o(l). 

By [6l Lemma 6] and the estimates of the first part of the paper, we have for any 
< r < 1/2 and for any p G n C, 

|||^°rGi,o(p)||e. <G||p||c, 
|||i?°r+'/'0ii„,+..(p)||6. s^GllplU.nc, 

\\\d\+'''g[^1_^M\\b.^g\\pU^^c- 

Splitting ID"! = |D0|^|D0|1-^ for some < t < 1/2, we thus get 

- Tr(|i?0|{Gi,o(«l^rcn), Gi;^o^___,(|(^n}) + Tr(|i^°|{Qiin, + ..(«;i^rcn), Gi,o(|^n}) - Oa^o(l). 

Finally, 

Tr(|i?"|{Qiin,+ («;z^rc„),ag2}) = --^7^Tr(|i?°|{Gi.o(«;z^rcn),Gi,o(|<^n})+Oa^o(«)- 

1 + aiJA 

We simplify the last expression in a general fashion. Let pi,p2 G n C and 
T(pi,p2) :=Tr(|i?0|{Gi,o(pi),Gi^o(p2)}). Then, 

nPi,P2) ^ 1^ f iE{p)+E{q))V;,ip-q)Vp,ip-q)Trc4M{p,q)M{q,p))dpdq 



1 

47r 



/B(0,A) 
S(0,2A) 



4./ i^AW^i^dfc 

is(0,2A) !«! 

BaD{pi, P2) - 4Tr f UA{\k\) 

JB(0.2A) 



Pl{k)p2{k) 

|fc|2 d'^- 



where we recall that 

„ . ^ ^_ [ 1 A _ l + (£ + fc/2).(^-fc/2) \ 

^ 7r2|fc|2 yK+fc/2KA £:(^ + fc/2) + £;(^-fc/2) ^ E{1 + k/2)E{£ - k/2) ) 

In our case pi = Kt'rcn and p2 = I'/'Pj since Uav € C because v Q CiC, we deduce 

Tr(|i?°|{Gi,o(«:i^rc„),Gi,o(|(/^n}) = SAi?(«:i^rcn, l"^!') + Oa^o(l) 
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thus 

1 + a±?A 

A.3. Estimate on Qs- As for Q2, we first write 
By the residuum formula, 

= \(p{KPren,-^)){S{R+)(p{KVren,-^-)\ + |S'(-R+)¥'(Kf ren, A)) ((^(KZ/ren, A) | , 

where ^ 

SO that \\S{R+)\\m^m < C||Q+(a)+Gi,o(Ki^re„)||c = Oa^o(l). Since ||Gg,_(g)||Q < 
CIIQIIq and Q+{a) - Q_(a) = + O^^oiV^), we have 

This leads to the expression 

Tr{\D°\{Qlin,+ {Kl/ren),aQ3}) = 2aRe((p(Kfren, A), { |D° | , Qlin,+ (k^ ren) }S'(-R+)</9(«;fc'ren, 
+ a TV (I £)0 I {Qlin,+ (K^^ren) , _ (| ^) (^| ) }) + O^^Q («) • 

Again we have 

|(<^(KZ/ren,A),{|D°|,Qii„,+ (Kt'ren)}5'(-R+)¥'(/«fc'ren, A))| = Oa^o(l), 

and 

Tr{\D°\{QnnA^u,en),Gi%_{\ip){cp\)}) = -Tr{\D°^^^^ 

+ T^{\D°\{Qii^,+A'^u,en),Gi%_^M{^\)})- 

Since ||giin,+,r(Ki^ren)||Q = Oa^o(l), ||G^1-,.(|V')(vI)IIq = Oa^o(l), (|<p) ||q 

O„^o(l), ||Gi,o(K!^ren)||e. = Oa^o(l), and \\\D^\Gi%_A\<p){<fi\)\\e, = Oa^o(l), 
we have 

-Tv{\D'>\{G,,oi^,.,,n),Gl,'^l_J\^}{^\)})+T^^^^ 

+ Tr(|D°|{Qnn, + ,r(«I^ren),G^l_^,(|¥')(¥'|)}) = Oa^o(l). 

A short computation shows that 

Tr(|7?°|{Gi,o(«I^ren),Go,l(|^)(¥'|)}) = - Tr(Gi,o («I^ren)i?|^) (^| ) , 

which implies 

|Tr(|i?0|{Gl,o(«I^ren),Go,l(k)(¥'|)})| ^ || Gl,0 (^^^^rcn) 1 1 6. 1 1 i?|^>(^| ! | e. 

||c||^||ii=Oa^0(l). 

Finally, 

Tr(|£>0|{Qii„,+ (/^iy,en),aQ3}) = Oa^oCa)- 
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The last term to estimate in (|2.68|) is thus ((^(Ki^rcn, A),{\D^\,Q-^r{a)}(p{Ki'rcn, A)). 
We do so by writing 

|((y5(KJ^ren, A),{|i:)°|,(3_,r(a)}</5(Kl'rcn, A))| ^ 2 1| (5_ (") 1 1 ffi |l 'y^('«^'rcn, A) || , 

and it remains to estimate ||(5-,r(Q:)|-D°|~^||i2_j.^2. RecaU that Q_.r(a) = Q-{a) — 
Q\in,-{KVrcn), which wc dcvelop up to 0{a). We have 

Q-,ria) = aGp'^i^ _(Gi,o('«t'rcn)) + aGi,o(p-(a) - BAUVrcn) + O^^gia), 
and writing for any Q that 

+^ I ^ — ^nAi?QnA-g — ^^nAK.,^„nA _ ^ — — 

271" Jr -D^ - fi^+iT] DX - Ai- + f / i'A - ^i-+^r| 

one obtains 

ll4?i,-(Q)l^°r'IU=^L= C||i?Q|i?°ri|U2^i2 = WRqWh^^l- ^ 2\\Q\\e, 
by the Hardy inequahty. Together with ||Gi,o(Ki/rcn)||62 — ^'a-i-oll); this imphes 

We also have ||Gi,o(p)||l2^l2 =^ ||Gi^o(p)||e2 ^ C'||p||c for any p, and thus 

l!Gi,o(P-(a) - SAKt'ron)|£'°r^||L2^L2 ^ || Gi,o(/0- (") " ^AKl^rcn) 1 1 L^^L^ 

< G||p_(a) - Bt,nv,,^\\c = Oa^o(l)- 

As a consequence, 

|((p(Kl're„,A),{|i:>°|,g_,^(a)}^(Kl'rcn,A))| = O„^o(a), 

and hence 

(A.2) Trop°Q+(a)) -Trop°Q-(a)) = (y^C^^zy^cn, A), i?V('*'^rc„, A)) 

~ ^y;(K;^i.cn, A), V;,^,,^„y)(Ki/i.cn, A)) +e'a_j.o(a). 

1 + Q;i?A 
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